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ABSTRACT 



Using stellar-dynamical calculations it is shown for the first time that counter-orbiting material emerges naturally in tidal interactions 
of disc galaxies. Model particles on both pro- and retrograde orbits can be formed as tidal debris in single encounters with disc 
galaxies of 1-to-l and 4-to-l mass ratios. A total of 74 model calculations are performed for a range of different initial parameters. 
Interactions include fly-by and merger cases. The fraction of counter-orbiting material produced varies over a wide range (from a 
few up to 50 percent). All fly-by models show a similar two-phase behaviour, with retrograde material forming first. Properties of the 
prograde and retrograde populations are extracted to make an observational discrimination possible. 

During such encounters the tidal debris occupies a certain region in phase space. In this material, tidal-dwarf galaxies may form. 
The modelling therefore can explain why galaxies may have dwarf galaxies orbiting counter to the bulk of their dwarf galaxies. An 
example is the Sculptor dwarf of the Milky Way, which orbits counter to the bulk of the disc of satellites. The modelling thus supports 
the scenario of the MW satellites being ancient tidal-dwarf galaxies formed from gaseous material stripped from another galaxy 
during an encounter with the young MW. 

A possible candidate for this galaxy is identified as the Magellanic Cloud progenitor galaxy. Its angular motion fits the angular motion 
of the MW disc of satellites objects. This scenario is in agreement with Lynden-Bell's original suggestion for the origin of the dSph 
satellites and the near-unbound orbit of the LMC. 

Key words, galaxies: kinematics and dynamics - galaxies: interactions - galaxies: dwarf - galaxies: formation - Local Group 



1. Introduction 

The early-recognised anisotropic distribution of Milky Way 
(MW) satellite galaxies naturally lead to the propositi on that they 
are phase-space cor related tidal dwarf galaxies dLvnden-Belll 
Il976t lKunkellll979l) . With the advent of cold-dark-matter cos- 
mology the satellites were however viewed as the unmerged 
remnants of the bottom-up hierarchical structure formation pro- 
cess. In (standard) cold dark matter (CDM) cosmology, galax- 
ies form through accretion and merging of smaller systems. This 
suggests a simple origin for the dwarf spheroidal (dSph) satellite 
galaxies of the Milky Way: they are believed to be dark matter 
dominated subhaloes, with a luminous matter content, that are 
gravitationally trapped in the Milky Way halo and have n ot yet 
merged completely with their host dWhite & Reesl [l978). The 
details of this origin, in contrast, are still disputed and not even 
the nature of the dominant dark m atter component in dSphs is 
well defined (Gil more et al.l [2007h . However, Zwicky's (1956) 
proposition that new d warf galaxies form when galaxies inter- 
act and his suggestion dZwickvll 19371) that galaxie s contain dark 
matte r today lead to the Fritz Zwicky Paradox ( Kroup a~et alj 
2010) since the implied standard cosmological model leads to 
a large population of TDGs that ought to have the properties of 
dE and dSph galaxies dOkazaki & Tanig uchi 2000). This would 
leave little room for the existence of cold-dark-matter dominated 
satellite dwarf galaxies, clashing with the expected existence of 
large numbers of such objects. 

In addition to this, the number of observed satellites was re- 
ported to be about an order of magnitude lower than expected 
from CDM simulations, resulting in what is termed the 'miss- 



ing satellite pr oblem' dMoore et al.l 119991: iKlvpin et all 119991: 
Diemande7aD|2008j). There are several proposed scenarios to 
solve this issue, like reducing the num ber of expected sate llites 
as only the most massiv e progenitors dLibeskind et al. l l2005l) . the 
earliest formed haloes ( Striga ri et al . 2007) or even only the low- 
mass systems dSales et d]l2007h~ are to be expected to produce 
satellite galaxies surviving until today. Other possibilities are to 
assu me that only the mo st massive subhaloes were able to form 
stars dStoehr et al.l 12002). All of these proposals are rather 'ad 
hoc'. The alternative proposed solution is that very faint satel- 
lites, that were until now over looked because of their low stellar 
densities (iTollerud et alj|2008l) . are still to be discovered. None 
of these proposals are satisfactory. 

Still, not only the observed number of Milky Way satellite 
galaxies is in conflict with the expectation, they also do not 
match the dark-matter mass - luminosity relation n or the form of 
the m ass function of luminous dark matter halos (Kroup a~et al.l 
20101) . The distribution is exceed i ngly anisotropi c, too (e.g. 



Lvnd en-Bellll 1 9761: iMaiewskill 1 994t [Hartwick 2000). Especially 



the most luminous 'classical' satellites pose a problem for hier- 
archical structure formation. They make up a significantly pro- 
nounced disc of satellites (DoS) , a structure highly inclined to 
the pl ane of the Milky Way disc dKroupa et alJl2005l:lMetz et all 
2007). The inclusion of more recentl y found, very fain t satellite 
galaxies supports this spatial feature (Metz et alj|2009l) . In fact, 
treating the 13 new ultra-faint satellites independently leads to 
virtually the sa me DoS solution as gi ven by the 1 1 bright classi- 
cal satellites (Kroupa et al. 2010). As Kroupa et alj d2010l) show, 
the standard CDM model (e.g. lLibeskind et al. 2009) can at best 
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Fig. 1. Orbital poles for the eight MW satellite galaxies for 
which proper motions are measured. The plot is an Aitoff pro- 
jection of galactocentric coordinates. The green lines give the 
uncertainties in direction of the orbital poles. The red star marks 
the mean direction of the orbital poles of the six satellites con- 
centrated in the centre, excluding Sagittarius and Sculptor, the 
red loop marks the spherical standard deviation. The dashed cir- 
cles are the regions 15° and 30° degrees from the normal of the 
DoS fitted to the 1 1 clas sical MW satellites. The figure is adapted 
from lMetz et ail (120081) . See the on-line edition of the article for 
a colour version of this figure. 



reproduce 0.4 percent of all MW-type dark matter halos having 
a host galaxy of similar luminosity as the MW with 1 1 satel- 
lites in a DoS. An i nfall of the satellite s as a group of dwarf 
galaxi es suggested bv lLi & Helmil (2008) and D'Onghia & Lake! 
J2009t) has been ruled out bv lMetz et a l. (2009) because of the 
thin structure of the DoS in comparison to the extension of ob- 
served dwarf galaxy groups. 

Investigatin g the available proper motions of the satellite 
galax ies (e.g.: iLvnden-Bell & Lvnden-Beill 119951: iPalma et"aT] 
2002, who also include globular clusters in their studies), the 
situation becomes even mo re puzzling . Out of eight satellite 
proper motions analysed by Metz et alj (l2008l) . six are found to 
be in agreement with being co-orbiting within the DoS, as vis- 
ible in their figure 1 (reproduced in Fig. Q] here). The authors 
found that such a situation is extremely unlikely even if the satel- 
lites were drawn from CDM simulations specifically tailored to 
solve the Milky Way DoS problem. Only two satellites have or- 
bital poles not directly associated with other poles. Of these, 
one (Sagittarius) is orbiting perpendicularly to the DoS, while 
Sculptor is on a counter-rotating orbit within the DoS. 

An alternative origin of the Milky Way satellite galaxies 
would be them forming as tidal dwarf galaxies (TDGs) perhaps 
about 10 Gyr ago in an interact i on of the young MW with an- 
other galaxy. ISawa & Fuiimotol d2005) proposed that the MW 
had an off-centre interaction with the primordial Andromeda 
galaxy M31 about 10 Gyr ago, whereby the M31 and MW 
satellite population was born. Another possibility would be a 
merger of two progenitor galaxies forming part of today's MW. 
The tidal expulsion of material during such an interaction might 
have led, through gravitational instabilities followed by clump- 
ing of matter within tidal tails , to the c reatio n of TDGs. This 
concept was first proposed by IZwickvl d 19561) and it was firs t 
observed to happen in the Antennae by Mira bel et alj ([1992). 
TDGs were subseque ntly found in many obs e rvations of inter- 
acting galaxies (e.g .: Hunsberger et alj [19961 : IWeilbache r~et~aT1 
12001 IWalter et al] 12006^ as well as in N -body calculations 
dWetzstein et al.l2007l: IBournaud & Due 2006). They woul d nat- 
urally lie in the plane of the interaction. lLvnden-Belll (fl976) first 
suggested that the Milky Way satellites might be of tidal origin. 



The common direction of orbit around the Milky Way, found 
for most satellite galaxies with measured proper motions, can be 
considered to be an indication for this scenario. The only excep- 
tions are Sagittarius and Sculptor. While Sagittarius is near to 
the Milky Way disc and thus may have suffered p recession or 
even a scattering encounter with another satellite (IZhaolll998l) 
and consequently can be expected to be far from the direction 
of the initial orbit, the Sculptor dwarf spheroidal galaxy is an- 
other case. The direction of Sculptor's orbitel angular momen- 
tum vector indicates that it indeed orbits within the DoS, but on a 
counter-rotating orbit. While this might, on the first glance, seem 
to pose a problem for the common origin of the DoS galaxies as 
TDGs, the opposing direction could in fact be another indication 
for this origin. If tidal material in galaxy encounters were to have 
only one orbital direction in most interactions, this would make 
a tidal origin of the MW satellite system unlikely. 

The aim of this contribution is to investigate whether 
counter-orbiting material can be formed in tidal interactions. 
Not the creation of actual TDGs is studied, but only the mo- 
tions of particles in N-Body models which represent stellar 
streams. Modelling of galaxy interactions and studies of the 
formation of TDGs in high-resolution calculati ons of gas-rich 
galaxy-encounters have prev i ously been don e (Bourn aud & Dud 
2006; IWetzstein et"all l2007t iBournaud et~ai1l2008bl) . But until 
now they have not focused on the orbital angular momenta of the 
created dwarf galaxies and the authors of this contribution know 
of no published work examining the possible creation of counter- 
rotating orbits. In the following, angular momenta of particles in 
interacting galaxies are analysed for the first time, searching for 
the creation of counter-orbiting material. 

In total 96 N-Body calculations including g as dynamics 
have been performed by IBournaud & Dud d2006l) . They iden- 
tified almost 600 substructures, but due to computing-time re- 
quirements the resolution was only sufficient to resolve objects 
more massive than 10 8 Mq, while all the satellite Galaxies of 
the MW except the Large and Small Magellanic Clouds have 
stella r masses in the range of 10 4 to 10 7 Mq (Striga rTet alj 
120081) . IBournaud & Dud d2006l) state that only a small frac- 
tion of their TDG candidates survive, a statement which has 
to be understood in the light of their detection criterion: their 
TDG candidates have to have a mass of at least 10 8 Mq. TDGs 
that loose mass and drop below this threshold are counted as 
non-surviving, whereas TDGs still exist (and can be formed) 
at much lower masses. Indeed, ex plicit high-reso lution calcula- 
tions of early dRecchi et al.1 12007 ) and long-term (Kroupa 1997; 
iKlessen & Kroupal 1 1998b evolution of 10 7 - 10 8 M TDGs 
demonstrate that th ey readily survive . Thus t he mass and spatial- 
resolution of the Bour naud & Du e (2006) models influences 
the deduced survival of TDGs and also their formation. In the 
highest-resolution calculation to date that includes gas, more 
than 100 young stellar objects, star clusters and TDGs, with 
masses between 10 5 an d 10 9 M(T), have been i dentified in one 
single galaxy collision dBournaud et alj|200 8bfl Furthermore, 
TDGs might also form from s tar cluster complexes wh ich are ob- 
served to exist in t idal ta ils dFellhauer & Kjoupal2002l) . Finally, 
Bournaud & Due (2006) performed merger-only calculations, 
whereas this contribution also includes fly-by interactions. This 
is why their results about the dependence of TDG survival- 



1 The more massive super star clusters identified in their work, with 
masses below 10 s Mq, can be interpreted to be co mpact spheroidal 
galaxie s in contrast to massive and rotating TDGs iBournaud et alj 
20081). 
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chance on the place of formation along the tidal tails can not 
simply be adopted herfl 

A full analysis of the orbital motions of formed TDGs would 
therefore demand calculations with very high resolutions that in- 
clude gas dynamics. Performing a parameter scan with these is 
currently not feasible. Gas dynamics is therefore neglected in 
the present work, to allow for a higher number of calculations 
scanning over initial velocities and directions between the inter- 
acting galaxies. This results in not looking for the formation of 
TDGs, but only in a study of the phase-space properties of the 
formed tails, testing if the concept of counter-orbiting material 
originating in a single encounter is feasible. 

This approach is applicable as angular momentum conser- 
vation implies that a condensed object (a TDG) will have the 
composite orbital angular momentum of the material (the parti- 
cles) accreted into it. While the effect of dynamical friction is 
stronger for a dwarf-galaxy- sized object th an for a single star, it 
is still negligi ble for TDGs (Kroupa 1 19971) . According to equa- 
tion 7-26 in dBinnev & Tremainel 1987b . the dynamical friction 
timescale, tfn C , for a M = 10 8 Mq galaxy orbiting at a galacto- 
centric distance d = 60kpc with a velocity of v = 125 kms~' is 
of the order of 10 2 Gyr, much longer than a Hubble time. Note 
that ff r i c oc thus the parameters chosen here are conservative, 
as the velocity v is low, d is small and fee will be even longer for 
TDGs forming with a variety of masses M < 10 8 M©. Therefore, 
the dynamics of TDGs and the stellar streams will not separate 
and the results can be interpreted in the light of a tidal origin for 
the MW satellites. 

Using a new versi on of the original (lMadeiskv&Bienll993t 
Fellh auer et al.l l200"ol) Superbox code developed by M. Metz, 
Superbox++, it is shown that particles on both pro- and retro- 
grade orbits are indeed a natural outcome of an encounter of 
two disc galaxies. Calculations are performed for perpendicu- 
larly oriented galaxies of 1:1 and 1:4 mass-ratios, scanning a 
range of different initial parameters and interactions including 
fly-by and merger cases. In almost all cases counter-orbiting ma- 
terial is found. 

Section|2]defines the nomenclature used in the following. In 
Sect. 13.11 the galaxy model, initial conditions and the thoughts 
leading to it are introduced. The results of two exemplary com- 
putations are described in Sect. [4] Section [5] gives an overview 
of the results from the other models calculated, scanning a range 
of initial parameters. A discussion including extraction of possi- 
ble observational constraints is performed in Sect. [6] where also 
the possibility of a tidal origin of the MW satellites system is dis- 
cussed. Note that Sect.|4]concentrates on equal-mass encounters, 
while Sects. [5]|6]generalize the encounters to non-equal galaxies. 
Finally there are concluding remarks in Sect. [7] There is an ap- 
pendix which gives details on how the scaling of the MW was 
performed. It also compiles the results of the parameter scans in 
12 tables. 



2. Definitions 

The galaxies in the models will be referred to as: 

Target galaxy: around this disc galaxy the particles ripped out 
of the infalling galaxy orbit and are later analysed. 



2 They find that TDGs formed at the tip of the tidal tail have the best 
chances to survive. However, in contrast to a merger where the tip will 
remain relatively undisturbed for a long time, that tip in a fly-by hits the 
second (target) galaxy first, which might lead to different conclusions. 



Infalling galaxy: the disc galaxy approaching the target. This 
galaxy will develop tidal tails from which particles will 
stream to the target and form retrograde and prograde or- 
biting material. 

Thus, in the following the terms prograde and retrograde will be 
used frequently: 

Prograde particles: particles that, after the galaxy-interaction, 
have an orbital angular momentum vector oriented in the 
same direction as the galaxies' orbital one within a circle 
of acceptance (COA) of a given angular radius. 

Retrograde particles: particles that, after the galaxy- 
interaction, have an orbital angular momentum vector 
oriented in the opposite direction than the galaxies' orbital 
one within a COA. 

Furthermore, proper motions of satellite galaxies in the vicin- 
ity of the MW are oriented preferentially in one direction, with 
the Sculptor dwarf galaxy moving along an opposite orbit. It is 
said that Sculptor is a 'retrograde' satellite, with respect to the 
bulk of the satellites' motions. This statement does not include 
any knowledge of a possible interaction with another galaxy and 
thus might cause some confusion. What is seen today as a 'ret- 
rograde dwarf satellite' might have been the prograde orienta- 
tion of the orbit when compared to the angular momenta of the 
galaxies in the initial interaction. Due to this ambiguity the use 
of 'prograde' and 'retrograde' is refrained from for those situa- 
tions where the initial interaction geometry is unknown, as is the 
case in the MW-satellite system. Instead we define: 

Co-orbiting: orbiting in the same direction as the bulk of the 
material (particles in the models or satellite galaxies in the 
MW). 

Counter-orbiting: orbiting in the opposite sense than the bulk 
motion. 

Thus, from the point of view of an observer today, Sculptor is 
counter-orbiting with respect to the other satellites. Finally the 
naming scheme of the models is defined: 

Names of models: The names of the models presented in the 
next sections are composed of the direction of the initial 
velocity vector for the galaxies, expressed in degrees from 
the z-axis (e.g. '7.5deg' meaning 7.5°) and the initial rela- 
tive velocity in per cent of the parabolic velocity for point- 
masses of the same mass (e.g. 'lOOvel' meaning 100 percent 

Of Vparab)- 



3. Models 

3. 1 . Scenario and galaxy model 

The models are based on the following scenario. About 10 Gyr 
ago, the progenitor of the MW had an interaction with another 
progenitor disc galaxy. During this interaction, tidal debris in 
the form of gas and possibly stars was formed and possibly the 
initial disc was destroyed in the process that might have included 
a merger. A merger of the two initial discs possibly formed part 
of the bulge of today's MW, whereas the MW disc seen today 

formed later through ongoing accretion. 

Investigating its chemical evolution, Ball ero et al.l (l2007h 
found that the bulge of the MW formed rapidly, with a timescale 
of the order of 0.1 Gyr. This is in agreement w ith a merger- 
induced origin. However, Babusiaux et al. (2010) suggests that 



3 



Pawlowski, Kroupa & de Boer: Making Counter-Orbiting Tidal Debris 
Table 1. Parameters of the dwarf disc galaxy models. 





equal mass ratio 


4-to- 1 mass ratio 




Symbol 


Infalling & Target 


Infalling 


Target 


Description 


Disc 2V d 


500000 


500 000 


500 000 


number of disc particles 


Mi 


8.0 x 10 9 M© 


4.0 x 10 9 M© 


16.0 X 10 9 Mq 


disc mass 


Ri 


1.60 kpc 


1.15 kpc 


2.25 kpc 


disc scale length 




8.00 kpc 


5.75 kpc 


11.25 kpc 


maximum disc radius 


Rq 


6.00 kpc 


4.30 kpc 


8.40 kpc 


solar radius 


Qo 


1.70 


1.70 


1.70 


Toomre Q at Rq 




0.25 kpc 


0.20 kpc 


0.35 kpc 


disc scale height 




2.50 kpc 


2.00 kpc 


3.50 kpc 


maximum particle height 


Halo A'h 


1000000 


1000000 


1000000 


number of halo particles 


M h 


8.0 x 10 10 Mq 


4.0 x 10 10 M© 


16.0 x 10 10 M© 


halo mass 


a 


10.0 kpc 


6.9 kpc 


12.4 kpc 


halo core radius 


nmiix 
K H 


100.0 kpc 


100.0 kpc 


150.0 kpc 


maximum halo radius 



the MW bulge consists of two parts: a spheroidal that might orig- 
inate from past merger events and a bar-like structure or pseu- 
dobulge. Thus, both a merger scenario forming a bulge as well 
as a fly-by interaction with less severe consequences for the early 
MW, but leading to the formation of a bar instability in the disc, 
appear plausible and are investigated in this contribution. 

To form a realistic bulge in the calculations would require the 
inclusion of gas, out of which the bulge stars can form within a 
short timescale. Furthermore, dissipative processes in the gas of 
galactic discs would lead to a change in the remnant-galaxy dy- 
namics compared to the purely stellar case. One therefore does 
not learn much from a comparison of the resulting central ob- 
ject in the calculations with the MW bulge, as the model only 
represents a solution to the collisionless Boltzmann Equation. 
Similarly, after the interaction the model disc is not comparable 
with today's MW disc as it, in this scenario, forms through ongo- 
ing accretion of intergalactic gas after the interaction took place. 
The present contribution therefore concentrates on an analysis 
of the general phase-space properties of the tidal material and 
does not claim to be a complete model for the formation of the 
MW disc, bulge and satellite system. Some remarks on the mor- 
phology of the remnant galaxy in the different interactions can 
be found in Section 

Assuming the MW bulge formed in an equal-mass merger 
gives approximately 0.5 MW bulge masses as a natural starting 
point for the mass of the early MW modelfl The disc galaxy pa- 
rameters are then scaled to that mass, while the dark matter halo 
parameters are adjusted such as to model the empirically found 
flat rotati on curves, the ori gin of which is a matter of ongoing 
debate dKroupa et al.ll2Q10D . The resulting parameters are com- 
piled in TableQ]for the equal-mass and 4-to-l mass models. The 
models for the 4-to-l mass ratio interactions were determined 
with the same method, but this time the mass was either half or 
twice the original value for target and infalling galaxy, respec- 
tively. More details on how these scalings were performed can 
be found in Appendix lAl 

Particle realisations of the models with 5 x 10 5 particles for 
the disc and 1 for the halo co mponent are set up using the 
code MaGalie (Boily et al. 2001). After that, the galaxy models 
are virialised by integrating them in isolation for 1 .5 Gyr. In the 



3 This might on the one hand over-estimate the mass of the classical 
bulge component of the MW in merger-interactions and on the other 
hand lead to a too-low pseudobulge-mass in case of fly-by interactions, 
but is considered to be a good compromise. 



case of the equal-mass interactions, after virialisation the model 
galaxy is duplicated to form the target and the infalling galaxy. 

One drawback of this approach is that it does not include 
the mass accretion forming the MW. The integrations cover 10 
Gyr and during this time the newly accreted disc should increase 
its mass until its present-day value. A changing mass will in- 
evitably change the orbital parameters of the pro- and retrograde 
particles, like the eccentricity and apocentre distance. But as the 
main aim is to show that counter-rotating orbits of tidal debris 
are in principle possible, that is no hindrance. The direction of 
the angular momenta of the particles will not be affected by a 
slow change, like that of a slowly forming disc, and their orbits 
will adjust adiabatically. 

3.2. Geometry and initial conditions 

This section focuses on describing the initial parameters of two 
exemplary equal-mass models of interactions between the two 
galaxies: one collision ending with the merging of the two disc 
galaxies (model 7.5degl00vel) and one off-centre collision lead- 
ing to a fly-by (model 5deg200vel). In both models, the infalling 
disc spins prograde with respect to the orbit of the two galaxies, 
because long lived tidal dwarfs form more easily in this geom- 
etry, as iBournaud & Dud d2006l) have stated. Results of further 
models scanning a range of initial parameters, including retro- 
grade spins of the infalling galaxy in the case of mergers, are 
summarised in Sect. In the case of retrograde fly-bys, tidal 
tails transferring matter are usually not as pronounced as in a 
prograde orientation of the infalling disc. 

The initial distance d m \ of the galaxies is chosen to be two 
times their dark matter halo maximum radii, d m i — 2 x ,R™ ax = 
200 kpc. Initial relative velocities are chosen with respect to the 
parabolic velocity v para b = 87 km s _I calculated from the total 
galaxy masses M talget and Mj n f a nm £ , and their initial distance d m i, 
assuming the galaxies (including haloes) are point masses. The 
relative velocities of the galaxies in models 7.5degl00vel and 
5deg200vel are 1.0 ■ Vparab and 2.0 • v para b, respectively. 

The two galaxies are put to +100 (infalling) and -100 kpc 
(target) along the z-axis. In the case of a merger, each is given 
half of the relative interaction velocity. Thus the centre of mass 
is at the origin at all times. In the fly-by case the target galaxy 
is assigned a smaller velocity than the infalling galaxy to make 
sure the target and the orbiting particles around it stay within the 
model volume. 
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Table 2. Initial conditions of the calculations 



Model: 


7.5degl00vel 


5dej 


;200vel 


Type: 


merger 


fly-by 


Galaxy: 


Target 


intalling 


Target 


intalling 


Total mass [10 u, Mo] 


8.8 


8.8 


8.8 


8.8 


Position [kpc] 










X 


0.0 


0.0 


0.0 


0.0 


y 


0.0 


0.0 


0.0 


0.0 


z 


-100.0 


+ 100.0 


-100.0 


+ 100.0 


Velocity [km s~ x ] 










v x 


0.0 


0.0 


0.0 


0.0 


Vy 


+5.68 


-5.68 


+3.03 


-12.13 


Vz 


+43.13 


-43.13 


+34.67 


-138.67 


p 


7.5° 


7.5° 


5.0° 


5.0° 


Ldisc 


+z 


+x 


+z 


+x 


r min [kpc] 




6.1 




8.3 


Time 




1000 time-steps = 500 Myr 



The mass of the whole galaxy (disc and halo) is given in 10 10 Mq. 
Position states the position of the centre of mass of each galaxy in kpc 
with respect to the centre of the model volume. Velocity does the same 
for the initial velocities in km s , fi is the angle between Vm and the 
z-axis. Ldi sc gives the directions of the angular momenta of the discs 
parallel to the axes, for example, +z means the vector points along the 
z-axis in positive direction. Finally, r„„„ is the perigalactic distance be- 
tween the two galaxies at their first passage. 



The direction of the velocity vectors and the perigalactic dis- 
tances during the first approach, r m ,„, can be found in Table [2] 
Also given are the direction of the internal (spin) angular mo- 
mentum vectors Ldi sc of the discs. The infalling galaxy's disc 
plane is perpendicular to the target's plane, which is motivated 
by th e possible Local Volume source of this galaxy (Met z et all 
2009) and lies in the orbital plane of the encounter. The infalling 
galaxy is always seen edge on from the target, in other words, it 
approaches the target like a circular saw. The geometry of the in- 
teraction is sketched in Fig. [2] The whole interaction takes place 
in the y-z-plane of the volume, in which the retrograde and pro- 
grade particle orbits lie. 

3.3. SUPERBOX++ calculations 

All models in this work are integrat ed using a new version 
of the well known code Superbox dMadeiskv & Bienl 119931: 
iFellhauer et al.ll2000l) . Superbox++ is a new implementation of 
a particle-mesh scheme with a hierarchy of high-resolution sub- 
grids used in Superbox. Differences in the realisation of the al- 
gorithm make the new code much more efficient than the for- 
mer. Superbox had been written in the Fortran language with 
a particular focus on the minimisation of usage of random ac- 
cess memory (RAM), which is no longer an issue. Superbox++ 
is implemented in the modern C++ programming language us- 
ing object oriented programming techniques. The algorithm has 
been developed with a focus on the performance of the code, but 
at the same time keeping memory consumption at a low level. 
Superbox++ makes optimal use of modern multi-core processor 
technologies. Superbox conserves the angular momentum well, 
especially for the high resolution used in these models, a s has 
been shown in tests performed by Fellhau er et al.l d2000i) . The 
maximum AL/L for the same grid-resolution was determined to 
be less than 0.05 percent after 1000 time-steps. 

For the models calculated in this work, the number of grid 
cells per dimension is 64, and the grid sizes are 10, 30 and 400 




Fig. 2. Schematics of the geometry of the interaction and the ini- 
tial paths of the galaxies. The orientation of the axes is given 
at the right, the dashed lines symbolise the direction of motion 
of the galaxies. The vectors L give the direction of angular mo- 
menta of the disc spins, indicated by the small arrows on the 
disc edges (labelled v 10t ).The vectors v; n j sketch the direction of 
the initial velocity of the galaxies, the angle between it and the 
z-axis is f3. 

kpc for the inner, middle and outer grid of the disc component, 
respectively. These values are the same as those given in lA.4l for 
the virilisation. The models are calculated for 10 Gyr in 20000 
steps. To get the best statistics to detect even small percentages 
of retrograde particles, positions and velocities of all disc parti- 
cles of the infalling galaxy are analysed. Snapshots are produced 
in time-step intervals of 50, corresponding to one snapshot every 
25 Myr. 

4. Results 

At first, a general description of two exemplary calculations is 
given, followed by a detailed analysis of the pro- and retrograde 
particle populations. In Sect. |5]the exploration of the larger pa- 
rameter space with more models analysed in an analogous man- 
ner is discussed. 

4. 1 . General description 

4.1.1. Model 7.5deg100vel: merger 

In this model, the galaxies pass by each other after 1 .55 Gyr at 
a perigalactic distance of « 6.1 kpc. After another 0.4 Gyr, the 
galaxies start to fall back towards each other and merge at about 
2.5 Gyr after the start of the integration. The target galaxy's disc 
is destroyed in the collision, a spheroidal central remnant galaxy 
is formed. A few hundred of the target's particles are expelled to 
distances up to 50 kpc, but no well-defined large tidal tails form 
out of the target galaxy. Thus, the following deals with particles 
from the infalling galaxy. 

The prograde particles (as defined in Sect. |2]l form in great 
number because they move in the same direction as the infalling 
galaxy with respect to the target. Thus, when the infalling galaxy 
swings around and merges, many particles in the outer part will 
follow this direction. Furthermore, a short tidal tail between the 
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two galaxies transports particles from the infalling to the target 
galaxy soon after the first passage. The majority of these is pro- 
grade. 

A part of the particles becomes retrograde when the two 
galaxies finally merge almost head on. Even though the cen- 
tres pass each other with a preferred, prograde orientation, some 
particles in the outer region will pass the combined centre of 
mass on the opposing side, leading to counter-orbiting particles. 
Furthermore, a significant fraction originates in the tip of the 
tidal tail pointing away from the target galaxy. That tail forms 
after the violent passing of the two galaxies and falls back to- 
wards the target together with the infalling galaxyQ 

4.1.2. Model 5deg200vel: fly-by 

In this model the galaxies approach each other until they pass by 
after about 1 Gyr with a minimum distance of 8.3 kpc. At 3 Gyr 
the distance approaches 240 kpc and the galaxies can be con- 
sidered well separated. Two tidal tails form out of the infalling 
galaxy, one pointing towards the target galaxy, the other into the 
opposite direction. The latter one is of no further interest, the 
particles in it later fall back towards the remaining infalling (but 
now receding) galaxy. 

The target galaxy gets disturbed, its disc heats up, thickens 
and a bar instability forms. However, in general it remains a ro- 
tating disc galaxy. While a pronounced spheroidal component 
that can be interpreted as a classical bulge is not formed, bars 
in disc galaxies are a typical formation mechanism of pseudob- 
ulges. Due to the interaction geometry, the spin angular momen- 
tum vector of the target disc lies within the plane of the interac- 
tion, no tidal tails are formed from the target galaxy. 

Via the tidal tail between the two galaxies, particles of the 
infalling galaxy's disc become stripped from their original host 
and stream towards the target galaxy, around whose centre they 
orbit from now on. The particles that are stripped in this first 
phase follow retrograde orbits. The second phase begins once the 
tidal tail, dragged along with the drifting away infalling galaxy, 
sweeps over the target galaxy's centre at about 2.5 Gyr. The now 
infalling particles pass the centre on the opposite side, resulting 
in a counter-rotating orbit, making these particles the prograde 
ones. This evolution is shown in two representative snapshots in 
Fig. [3] Because they fall in from further away - the distance be- 
tween the galaxies has increased during phase one - it can qual- 
itatively be expected that prograde particles have higher orbital 
eccentricities and higher apocentre distances on average. 

4.2. Analysis: snapshots used, spherical coordinate system, 
choice of analysed particles 

It is the key issue of this work to analyse the movement of par- 
ticles around the target galaxy. As an analysis of the orbits of 
thousands of particles is not feasible, a simple and easy to visu- 
alise way to trace the direction of the orbiting particles is to look 
at their orbital angular momentum vectors L. Prograde and ret- 
rograde particles should show up as two separate populations in 
a plot of orbital angular momenta directions. That way the par- 
ticles can be identified, allowing to compare other properties of 
the two populations, finding differences that might be observa- 
tionally investigated. This makes the approach testable in nature. 



4 The two models are also visualised in movies, 
available online at the Downloads Page of the AlfA: 
http://www.astro.uni-bonn.de/download/animations/ 



The analysis of the orbits of the infalling galaxy's parti- 
cles caught by the target starts at 3 Gyr into the computation. 
Determinations of the angular momenta are carried out for every 
0.5 Gyr, until the end of the integrations at 10 Gyr. This results 
in 15 analysed steps per model integration. To determine the ec- 
centricities, every snapshot from the above mentioned start is 
used, resulting in a temporal resolution of the particles paths of 
25Myr. 

The orbital angular momenta are determined with respect to 
the centre of density of the target galaxy in case of fly-by inter- 
actions. For mergers, the equal-weighted mean of the infalling- 
and the target-galaxy centres of density is used. To display the 
the orbital angular momentum vectors L, a spherical coordinate 
system is set as follows: the angle 6* runs from (positive z-axis) 
to 180 degrees (negative z-axis). The angle <p gives the direction 
around this axis in a range from to 360 degrees. </> = 0° implies 
L has a negative y-axis component only and no x-component, 
whereas <p = 90° has no y- but solely a positive x-axis compo- 
nent. Thus, 6-0° points into the same direction as the spin 
angular momentum vector of the target galaxy's disc rotation. 
Particles exactly prograde with respect to the motion of the in- 
falling galaxy will have L pointing along the positive x-axis, thus 
to (6, (f>) = (90, 90). Consequently, retrograde particles have L 
pointing to (6, (p) = (90, 270). 

4.3. Orbital angular momenta, particle numbers and ratios 

The final distribution of orbital angular momenta distributions 
for the merger model 7.5degl00vel is plotted in the left panel 
of Fig. [4] For the sake of clarity only a subset of 20000 parti- 
cles is shown. The right panel plots the angular momenta for all 
particles within 400 kpc of the target galaxy in the fly-by model 
5deg200vel. The angular momenta in the merger model show a 
wider distribution than in the fly-by, but in both cases two pop- 
ulations near the expected positions of pro- and retrograde ma- 
terial can be clearly seen. Until the end of the calculations, the 
debris is distributed in a disc. This can be seen in Fig. [5] where 
the tidal debris is plotted in a projection edge-on to the orbital 
plane of the galaxies. 

To identify the pro- and retrograde particles, a circular area 
(circle of acceptance, COA) around the expected positions at 
(0, <p) pm = (90, 90) for a prograde and (0, 0) retro = (90, 270) for a 
retrograde orientation is looked at. The circular area has a radius 
of 30° in the case of a fly-by and of 60° in the case of a merger, to 
account for the stronger changes in angular momenta directions 
due to precession and torquing during the merger process in the 
latter case. 

To check whether most particles attracted by the target are 
either in the pro- or in the retrograde COA, the situation at the 
end of the calculations is looked at closer. For this, all particles 
with radial distances r < 400 kpc from the target at the end of 
the calculations (after 10 Gyr) are counted. For the merger model 
7.5degl00vel,2.1 xlO 5 pro- and 1.6xl0 5 retrograde particles are 
counted in the respective circles of acceptance. Relative to the 
total number of 5.0 x 10 5 particles considered, material outside 
the COAs accounts for merely 26 per cent. For the fly-by model 
5deg200vel, 4500 pro- and 3000 retrograde particles are counted 
from a total of 9100 particles within 400 kpc of the centre of 
density of the target galaxy. Thus, only 17 per cent are neither in 
the pro- nor in the retrograde COA. 

Figure|6]shows the number of particles as a function of radial 
distance from the target galaxy at the end of the integrations for 
both models. The relative amount of retrograde particles Pretro is 
expressed as the fraction of the number of particles N retro that are 
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Fly-by model 5deg200vel at 2.5 Gyr Fly-by model 5deg200vel at 3.5 Gyr 




Z in kpc Z in kpc 

Fig. 3. Time evolution of a projection onto the z-y-plane of the particle positions for model 5deg200vel, a typical fly-by. The particle 
numbers are logarithmically binned. Both plots are centred on the target galaxy. In the left panel the situation at 2.5 Gyr is shown, 
particles of the infalling galaxy (outside the upper left corner) stream towards the target galaxy (first phase). They pass the target's 
centre in the direction of more negative y-axis, defining their orbit to be retrograde. The right panel shows the situation in the second 
phase at 3.5 Gyr. Now the tidal tail has swept over the centre of the target galaxy and the particles still streaming towards the target 
pass it on the side of less negative y-axis. They end up on counter-rotating orbits compared to the particles of the first phase. Thus, 
particles on prograde orbits are formed after the retrograde ones in fly-by encounters. 
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Fig. 4. Directions of the orbital angular momenta of particles of the infalling galaxy for merger model 7.5degl00vel (left) and 
fly-by model 5deg200vel (right) at the final snapshot after 10 Gyr of calculation. Particles which have prograde motions have 
0° < <p < 180°, those with retrograde motions have 180° < (p < 360°. In both cases two distinct and significant populations, one pro- 
and one retrograde, show up. 8 = 0° points towards the target galaxy spin direction ("north"). Merger: the 20000 plotted particles 
were randomly chosen from all 10 5 particles with radial distances kpc < r < 170 kpc from the centre of the merged galaxies. 
They show a wide distribution, characteristic for the more chaotic motions after the merger. It can be seen that the bulk of pro- and 
retrograde particles have precessed away from the expected positions at 6 — 90°, <p = 90° or cp = 270° respectively, but still lie 
within the 60-degree circles-of-acceptance (COA) for mergers. Fly-by: the plot includes all particles which have radial distances of 
less than 400 kpc from the target galaxy. Both the particles with prograde and the ones with retrograde motion remain close to the 
expected orbital angular momenta. They lie well within the 30-degree COAs for fly-bys. Compare to Fig.[[2]which plots the MW 
satellite galaxy distribution in a similar coordinate system. 



found to be retrograde compared to the sum of pro- (N plo ) and 
retrograde particles, 



*Y pro N [- e 



(i) 



The fraction of particles that are outside of the COAs is given 
relative to the total number of particles 7V a n as 



Wall - (AW + A^retro) 



/COAs 



Wall 



It is highest near the centre for both models. Low apocentre par- 
ticles will later be excluded as they suffer most from interactions 
in the central part of the target galaxy. 

Particles with distances up to 160 kpc exist in the merger 
model. The majority of the infalling galaxy has merged with the 
target and is thus found at the centre. P letro varies with radius, 
starting with an equal distribution at low radial distances, drop- 
ping below 0.25 for intermediate distances and rising again to 
a maximum of 1.0 at the highest distance of 160 kpc. For ra- 
dial distances larger than about 30 kpc /qq^ is less than 10 per 
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Fig. 5. Edge-on views of the tidal debris discs in the merger model 7.5degl00vel (left panel) and the fly-by model 5deg200vel (right 
panel) at the end of the calculations. Retrograde particles are plotted in red onto prograde particles (green), both are plotted on 
top of 20000 particles arbitrarily selected from the target galaxy (blue). In the merger case, only the first 20000 particles of the 
pro- and retrograde population are plotted. The tidal debris form a disc around the central galaxy which is oriented in the plane 
of the interaction. Note that the thickness of the tidal debris distribution is less than the thickness of the MW DoS (Met z et al.l 
120071: iRro upa et al.l 1201 Oh . In a less-ideal interaction geometry (e.g. when the infalling disc or the galaxy-orbit is not perfectly 
perpendicular to the target disc) and through subsequent dynamical evolution, the debris-disc can be expected to be thicker. See the 
on-line edition of the article for a colour version of this figure. 



cent. It rises again at the outermost bins, where the total particle 
numbers are of the order of tens only and velocities are slowest. 

In contrast, in the fly-by model 5deg200vel, most particles 
are within 130 kpc, but there is a constant value of approximately 
30^40 particles per 16-kpc-bin at larger distances. Particles 
bridge the whole gap between target and infalling galaxy, which 
have a separation of 720 kpc by the end of the calculation. The 
transition from retrograde particles produced in phase one of the 
encounter to prograde particles from phase two can be clearly 
seen. Outside of 80 kpc, no retrograde particles are found. P re tro 
consequently drops to zero, but has a maximum of 0.75 in the in- 
ner 50 kpc. As in the merger case, the vast majority of particles 
is inside the pro- or retrograde COA and /jSoAs e is again highest 
in the central region, but still reasonably low at less than 25 per 
cent. 

To be accepted for the detailed analysis, particles have to 
approach both an apo- and a pericentre within the analysis-time 
between 3 and 10 Gyr. Those particles which are found within 
the respective (pro- and retrograde) COA in at least 12 out of 
the 15 analysed snapshots are then identified as prograde and 
retrograde particleo Finally, the apocentre distance has to be at 
least 30 kpc. Thus the numbers of included particles are lower 
limits only. The benefits, however, outweight this: all accepted 
particles are on bound orbits around the target galaxy and can be 
assigned an orbital eccentricity (see Sect. |4.4| for details). 

5 The direction of L might change in the course of time, especially 
particles of highly eccentric trajectories can be strongly perturbed by lo- 
cal gradients in the potential. Furthermore, as the potential is not spher- 
ically symmetric, the path of the particles might be changed slightly 
on their orbit by precession and torques, changing their orbital angular 
momentum vectors and moving them out of the COA. Thus, demanding 
the particles to appear inside the COA for all snapshots was deemed to 
be unrealistic. 



In the case of merger model 7.5degl00vel, there are 36800 
prograde and 6000 retrograde particles with apocentre distances 
of at least 30 kpc. The fraction of retrograde particles with re- 
spect to the sum of pro- and retrograde ones is P retr o = 0.14. 
Particles with apocentres larger than 60 kpc are considered to be 
'far'. For these, the numbers are: 15000 prograde and 2250 ret- 
rograde particles, the fraction of retrograde particles with apoc- 
entres above 60 kpc is P{^ tIO = 0.13. 

For fly-by model 5deg200vel, the numbers are 3450 pro- 
grade and 2300 retrograde particles with apocentre distances 
> 30 kpc. In this case P re ti- is 0.40. The numbers for 'far' par- 
ticles in this case are: 1920 pro- and 170 retrograde, leading to 
a lower retrograde fraction P f r f tIO - 0.08, as expected from the 
two-phase origin in the fly-by. 

In the two models, significant numbers of particles orbiting 
in both directions are found. This is, as far as the authors of the 
present contribution know, the first time such an analysis was 
performed, showing that it is indeed possible to form particles 
on opposed orbits in both, mergers and close passings of two 
galaxies. This shows that, as tidal dwarf satellites form from 
such material and thus share the same phase-space region, co- 
and counter-rotating orbits of TDGs can in principle be created 
in a single galactic encounter. 

4.4. Eccentricities and apocentre distances 

The particles identified a posteriori as pro- and retrograde ones 
are now analysed in detail. For this analysis, all particles that are 
found within the afore mentioned 60- or 30-degree circles of ac- 
ceptance in at least 12 analysed steps are taken into account. This 
number of steps results in particle numbers closest to the ones 
derived from averaging over the results with minimum numbers 
of 10 to 15 demanded steps. 



8 



Pawlowski, Kroupa & de Boer: Making Counter-Orbiting Tidal Debris 



10° 



Merger model 7.5degl00vel 



10 u 



- - all 

— prograde 
retrograde 




^ 10 



20 40 60 80 100 120 140 160 180 




Fly-by model 5deg200vel 



- - all 

— prograde 
retrograde 



50 100 150 200 250 300 350 400 




80 100 120 140 160 180 

r in kpc 




100 150 200 250 300 350 400 




150 200 250 300 350 400 

r in kpc 



Fig. 6. The radial distribution of particles in the merger model 7.5degl00vel (left) and fly-by model 5deg200vel (right) at the final 
snapshot. The top panels give the number of particles in each bin, the 25 bins are plotted logarithmically to account for the vast 
range in numbers. They have a width of 6.8 kpc in the merger and 16 kpc in the fly-by case. Particles counted as prograde and 
retrograde are those within the respective COAs of 60° (merger) and 30° (fly-by) radius. The middle panel shows Pretro, the fraction 
of retrograde material. The bottom panel plots J?^f e , the fraction of all particles that are neither in the pro- nor in the retrograde 
COA. Merger: both pro- and retrograde material exists out to a maximum distance of about 160 kpc. P re tro is highest in the centre 
and at the maximum distance, varying between approximately 20 and 40 per cent. f^ff^T is highest in the central region, where 
the most chaotic motions are expected, but particles outside of the COAs are a minority at all distances. Fly-by: while retrograde 
particles are strongly centred and do not exist beyond about 80 kpc, particles on prograde orbits can be found all the way between 
the two galaxies. Thus, a fly-by can in principle produce prograde tidal debris at all distances between the interacting galaxies. The 
plot of Pretio shows that only in the inner region retrograde material dominates. As it is formed early, when the two galaxies are 
still closer to each other, it has lower average apocentres. The lower panel of /qq^ shows that the majority of all material is found 
within the COAs, as supported by the strong concentration of the particles in Fig. |4] (right panel). The maximum is found in the 
inner region, where about 25 per cent of the particles are outside the COAs. See the on-line edition of the article for a colour version 
of this figure. 



To describe the elongation of particle orbits each particle is 
assigned an orbit 'eccentricity' e by determining the apocentre- 
4spo, and pericentre distances, af per i, of its orbit, 

^apo — ^peri 

e — . 

^apo "t" ^peri 

To determine c/ apo and d per i the particle trajectories, in a coordi- 
nate system centred either on the target galaxy's centre of density 
(for fly-bys) or on the average of the centres of density of target 
and infalling galaxy (for mergers), is followed starting with the 
last snapshot and proceeding backwards in time. The pericentre 



distance d pel -i is estimated from a circle fitted to the three time- 
steps (±25 Myr) closest to the minimum distance of the particle. 
This avoids an overestimation of d pa -\ which would result in an 
underestimation of e in case of fast particles on highly eccentric 
orbits which have a coarse trajectory sampling. 

In the case of the merger model 7.5degl00vel the orbit ec- 
centricities stretch from low values, 0.5 for pro-, 0.35 for retro- 
grade particles, up to an eccentricity near 1 .0, which describes a 
parabolic orbit. The prograde particles are centred around their 
mean at e pro = 0.76. The standard deviation of the distribution 
has a value of 0.09. The orbit eccentricities of the retrograde par- 
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Fig. 7. Apocentre distances for merger model 7.5degl00vel (left) and fly-by model 5deg200vel (right). Plotted are the numbers of 
particles in bins of <4 po with width 4.0 kpc, starting at 30 kpc. In case of a merger, particles on pro- and retrograde orbits (within 
the respective COAs ) can be found at all apocentre distances from 30 up to about 140 kpc. In the fly-by case, in contrast, pro- and 
retrograde particles have distinct apocentre properties. The retrograde particles have lower apocentres (up to about 80 kpc) while the 
prograde particles have maximum apocentre distances of 120 kpc, but only few of them are found with d. dpo below 50 kpc. This hints 
at the two-phase origin typical for fly-by interactions. Particles with retrograde orbits are formed first, while the passing galaxies 
are still close to each other, and thus have lower average apocentre distances. 



tides in contrast are distributed more widely. Their mean value 
is e re tro = 0.76 with a standard deviation of 0.18. Both particle 
groups have apocentre distributions, plotted in Fig. [7] populated 
mostly near the minimum accepted distance of 30 kpc, beyond 
which their numbers drop. Particles of both groups reach out to 
140 kpc. There thus is co- and counter-orbiting material at sim- 
ilar radial distances from the target galaxy. The average eccen- 
tricity of the subset of 'far' prograde particles (t/ apo > 60 kpc) is 
essentially the same as for 'all' prograde particles (0.77), but it 
is higher (0.94) for 'far' retrograde particles. 



The situation is different for the fly-by model 5deg200vel. 
As expected, the later-infalling prograde particles have high ec- 
centricities, with a minimum value of about 0.65. Their mean 
value is <? pro = 0.82, around which the distribution is concen- 
trated, the standard deviation is only 0.08. The retrograde parti- 
cles have a minimum orbit eccentricity of about 0.5, the mean 
is e re tro = 0.80 and the standard deviation is 0.11. As in the 
merger case, with 0.93, the average eccentricity of the subset of 
'far' retrograde particles is higher, while there is no significant 
change in average eccentricity when looking only at 'far' pro- 
grade particles. The apocentre distance distribution of the fly-by 
model is shown in the right panel of Fig. Q In this plot, prograde 
and retrograde particles have clearly distinguishable properties. 
The prograde orbits have apocentres in the range between the 
accepted minimum of 30 to approximately 120 kpc, with a max- 
imum near 60 kpc. Note that, as obvious from the right panel of 
Fig. [6] prograde particles spread out to over 400 kpc. They are 
not included in Fig. Q because they are still falling in towards 
the target galaxy at the end of the calculation, so for them no 
eccentricity is determined and they thus are excluded from the 
detailed analysis. In contrast to the prograde particles, the retro- 
grade particles are distributed more closely to the centre of the 
target galaxy. Retrograde particles of a significant number only 
appear up to about 80 kpc. 



4.5. Site of origin within the precursor disc 

For the particles found to be pro- and retrograde, the site of ori- 
gin within the infalling galaxy is determined. This is the average 
radial position measured from the centre of the infalling galaxy 
before it interacts with the target galaxy. In Fig.[8]histogrammes 
for both models show the distribution of initial distances for the 
pro- and retrograde particles respectively. 

In merger model 7.5degl00vel, particles have initial dis- 
tances as close as 3 kpc from the centre, the mean values D al1 
for both populations are similar: D'^ =5.9 kpc for the prograde 

particles and D^ - 6.0 kpc for the retrogrades. The widths 
of the distributions are similar, too. This shows that co- and 
counter-orbiting tidal derbis in this merger come from similar 
regions in the infalling galaxy. 

In fly-by model 5deg200vel the particles originate from a 
similar region in radial distance, spreading from about 4 kpc to 
the edge of the initial disc at around 8 kpcQ. However, the distri- 
bution within this region is different for pro- and retrograde par- 
ticles. The former originates from further inside the disc when 
compared to the retrograde particles, they show a peak at around 
4.5 kpc (and a second, lower but wider one at 6.5 kpc) in contrast 
to the retrograde particles' peak above 6.5 kpc. The mean values 
of the distributions are D^ Q =5.9 kpc in case of the prograde 

particles and D^ ti0 = 6.7 kpc for the retrograde ones. Thus the 
mean initial distances differ by about 800 pc. 

As visible in Fig. [3] the prograde particles appear later than 
the retrograde ones. They appear when the tidal tail of the in- 
falling galaxy sweeps over the centre of the target galaxy. The 
outer part of the tidal tail, near it's tip, consists of particles that 
are more weakly bound to the infalling galaxy and thus originate 
from the rim of the disc. This part of the tidal tail forms retro- 
grade particles. Correspondingly the prograde particles, formed 
further down the tidal tail, originate from closer to the infalling 
galaxy's disc centre. 

6 Some particles have an origin beyond the initial cut-off radius of 
8.0 kpc because the disc slightly expands during virialisation. 
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Fig. 8. Histogrammes of initial distances dim within the infalling galaxy for merger model 7.5degl00vel (left) and fly-by model 
5deg200vel (right) in bins of 0.25 kpc width. All particles within the respective COAs with apocentres above 30 kpc are included. 
In the merger case, the particles of both populations initially had similar radial distances from the centre of the infalling disc galaxy. 
In the fly-by case, the particles on retrograde orbits had on average higher initial distances from the infalling disc galaxy's centre. 
Thus, they originate more closer to the rim of the galaxy. There are particles on prograde orbits coming from a similar region, but in 
additon to those many come from further inside the infalling's disc. Thus, on average the prograde material originates from a more 
central part of the infalling galaxy. 



5. Exploration of the parameter space 

Varying the initial relative velocities and their directions between 
the two galaxies produces a wide range in all analysed proper- 
ties. For fly-by models, 16 calculations with four different ini- 
tial velocities (1.8, 2.0, 2.2 and 2.4 times the parabolic velocity 
v P arab) an d f° ur different initial velocity directions (5.0, 6.0, 7.0 
and 8.0 degrees from the line connecting the two galaxies, which 
is the z-axis) are performed. Mergers are set up with both pro- 
and retrograde oriented infalling galaxy discs, in total there are 
18 for each. They have initial velocities of 0.5, 1.0 and 1.5 times 
Vparab an d directions of (head-on collision), 2.5, 5.0, 7.5, 10.0, 
15.0 and 20.0 degrees from the z-axis. The fastest models only 
include directions up to 7.5°. Furthermore, several models with 
unequal galaxy masses are calculated. Their mass ratio is 4 : 1 
for Target : Infalling. 10 fly-by models and 6 merger models of 
each orientation are analysed. 

The values extracted from the calculations (e.g. particle 
numbers and retrograde-fraction, average initial distances, or- 
bit eccentricities) are compiled in Tables IB. 11 to IB. 121 in the 
Appendix for each model type (equal mass fly-by, both orien- 
tations in mergers and the same for the 4-to-l mass fraction in- 
teractions). 

Before discussing the results for the different types of inter- 
actions in detail, a general remark on the expected radial distri- 
butions of tidal material is asked for. A fly-by can, in principle, 
produce tidal debris at all distances, see Fig. [6] The tidal bridge 
between the two departing galaxies becomes stretched and will 
thin out, but the tidal debris near its centre has a small velocity 
dispersion and will feel both the attracting gravitational force of 
the infalling and of the target galaxy. Thus, it might stay near 
the centre of gravity of the interactions while the two precursor 
galaxies recede. A merger, in contrast, can only produce tidal 
debris up to distances determined by the energy of the merger 
process. 

Faster initial velocities lead to later mergers because the 
galaxies oscillate more strongly before settling to their common 
centre of mass. Therefore, the analysis starts later for models 



with high initial velocities. This delay is reflected in the max- 
imum orbital periods, they are shorter, reducing the maximum 
apocentre the analysis can find. 

Before discussing the different interaction types, some gen- 
eral remarks on the typical remnant galaxy shall be made. Its 
morphology depends on the type of interaction. In the case of 
equal-mass mergers the discs get destroyed and a spheroidal, 
bulge-like remnant forms. In contrast to this, 4-to-l mergers 
leave the more massive target galaxy disc identifiable, but it 
heats up and gets tilted. The inclusion of gas dissipative pro- 
cessed would probably help the formation of a spheroidal rem- 
nant. The target disc is preserved in fly-by interactions of both 
mass-ratios, however it is disturbed. It heats up and a bar insta- 
bility forms, which mig ht result in a pseudo bulge-component as 
observed in the MW (Ba busiaux et al.ll2010l) . 

The spin angular momentum vector of the target galaxy disc 
lies within the plane of the interaction in all models. Due to this 
interaction geometry, no pronounced tidal tails form out of the 
target galaxy. Therefore, the analysis concentrates on the parti- 
cles of the infalling galaxy. 

5. 1 . Equal-mass fly-bys 

Particles in both orbital directions are found in all analysed mod- 
els. The calculated fly-by models all show the same simple two- 
phase behaviour: once the two galaxies pass, a tidal tail evolves 
from the infalling galaxy. Falling towards the target galaxy, the 
tail first forms retrograde particles (phase one). After some time, 
as the galaxies depart further, the tidal tail sweeps over the centre 
of the target galaxy. This changes the orbital direction of the par- 
ticles falling onto the target afterwards, making them prograde 
(phase two). 

The total number of pro- and retrograde particles within the 
COAs for which apo- and pericentre were determined, iVaii = 
N pm + Metro lies m the range of about 800 to some 6000 particles. 
Up to 2000 particles reach apocentre distances above 60 kpc, the 
lowest number found there is about 400. 
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Table 3. Pro-and retrograde particle properties 



Equal-mass fly-by models 



Symbol 


all(d apo > 30 kpc) 


far (4 po > 60 kpc) 




prograde 


retrograde 


prograde retrograde 


7.5degl00vel 










(Merger) 










N 


36800 


6000 


15000 


2250 


f N [%] 


7.36 


1.20 


3.0 


0.45 


M [10 s M©] 


5.89 


0.96 


2.40 


0.36 


P 

1 retro 




0.14 


0.13 




e 


0.76 


0.76 


0.77 


0.94 


Ae 


0.09 


0.18 


0.08 


0.02 


< d ml > [kpc] 


5.89 


6.02 


6.75 


5.38 


A < dm > [kpc] 


1.33 


1.27 


0.95 


0.86 


rf™ d [kpc] 


5.96 


5.92 


6.85 


5.20 


5deg200vel 










(Fly-by) 










N 


3450 


2300 


1920 


170 


fs [%] 


0.69 


0.46 


0.384 


0.034 


M [10 7 M©] 


5.52 


3.68 


3.07 


0.27 


P 

1 retro 




0.40 


0.08 




e 


0.82 


0.80 


0.80 


0.93 


Ae 


0.08 


0.11 


0.08 


0.03 


< <4ni > [kpc] 


5.89 


6.73 


6.07 


7.23 


A < dj„i > [kpc] 


1.05 


0.68 


1.07 


0.71 


«C [kpc] 


6.02 


6.68 


6.23 


7.30 



N is the number of particles with measured apocentre of at least 30 (all) 
and 60 (far) kpc, corresponding to a total mass M and f N , the fraction of 
total particle number of the infalling galaxy. The mean orbit eccentricity 
is denoted with e, having a standard deviation of Ae. Similarly, < d ir[i > 
stands for the mean initial radial distance within the infalling galaxy, 
given in kpc. A < d m [ > is the standard deviation and rf™ cd the median 
value. The values were determined by demanding the particles to be 
inside the circles of acceptance (COA) for at least 12 of the 15 analysed 
steps. 



The fraction of retrograde particles, P re t ro (Eq.Q~Jl, can be seen 
in Fig. [9] It shows a clear trend: the faster the galaxies pass, 
the higher is P re tro- In a slow passage, prograde particles domi- 
nate, while in the fastest passages calculated, the retrograde par- 
ticles are the dominating population. This seems to be almost 
independent of r m ; n , the minimum distance of the two galaxies 
during their passage. The rmin-fretro-plane is quite evenly popu- 
lated, with the region of small r m i n and large P le tro being empty, 
^retro ranges from 0.09 in model 8degl80vel, the slowest pas- 
sage with high galaxy-galaxy-distance at pericentre, to 0.99 in 
model 8deg240vel, the model with the highest initial velocity 
and farthest pericentre distance between the galaxies. 

Considering only 'far' particles with apocentres above 60 
kpc (lower panel in the figure), this trend is still seen, but 
most fractions drop strongly. For the slowest passages (rela- 
tive velocity of 1.8 times the parabolic velocity), virtually no 
retrograde particles with higher apocentre distances are found. 
For the fastest models, in contrast, Plf lm only changes slightly. 
Therefore, P£ m still shows a wide spread, from 0.00 to 0.99. 



However, all but the fastest passage models have P[^ tm < 0.2. It 
thus seems that more extreme ratios of co- and counter-orbiting 
material are preferred in fly-by scenarios. The reason for the 
drop in P retro when considering only far particles for most mod- 
els is the sharp transition between the two phases, producing 
retrograde particles with low apocentre distances first. Many of 
these do not reach out to over 60 kpc. 
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Fig. 9. The fraction of particles on retrograde orbits after the fly- 
by interaction of two equal-mass galaxies. They have apocen- 
tres above 30 kpc (Pretio, upper panel), and above 60 kpc CPj^L,, 
lower panel). See Tables IbTI and HOI for the values. The frac- 
tion of retrograde material is well defined by the initial velocity 
only, without much dependence on the minimum distance rmin 
at which the two galaxies pass each other. P re tro increases with 
increasing initial velocity, which ranges from 1.8 ('180vel') to 
2.4 ('240vel') times the the parabolic velocity. See the on-line 
edition of the article for a colour version of this figure. 



It can be concluded that the relative velocity of the two in- 
teracting galaxies has a dominating effect on the ratio of pro- 
grade to retrograde particles. The perigalactic distance during 
the passage has a smaller effect. This is of importance for the re- 
construction of galaxy-encounter histories. Within the scenario 
of an early, similar-mass fly-by event, the fraction of retrograde 
material (e.g. in the form of satellite galaxies) gives direct hints 
at the relative velocity at which this encounter happened. But 
one has to keep in mind that there is an ambiguity: pro- and ret- 
rograde orbits in the models can not directly be identified with 
co- or counter-rotating orbits. This ambiguity has to be broken 
with additional information. For this, the use of the radial (apoc- 
entre) distribution of the fraction of counter-orbiting tidal debris 
is suggested. If it decreases strongly for larger distances from the 
central galaxy, the debris could be identified with the retrograde 
material in the models which shows the same behaviour. This 
information then would not only break the ambiguity, but also 
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Fig. 10. Prograde equal-mass fly-bys: the mean initial distances 
to the infalling galaxy's centre, and Df^ tIO , for particles on 
pro- and retrograde orbits after the interaction. Upper panel: all 
particles with apocentre distances above 30 kpc are taken into 
account. There is a clear trend that retrograde particles originate 
from further out. Lower panel: the mean initial distances Djjjjj, 

and £>fetro snow a better distinction, the before mentioned trend 
is more pronounced when only considering particles with apoc- 
entres above 60 kpc. All models of Tables IB. 21 were included, 
except models 6degl80vel, 7degl80vel and 8degl80vel as there 
are no far retrograde particles found. 



allow a reconstruction of the orientation of the orbital angular 
momentum of the early encounter. 

Plotting the average initial distance D versus pericentre pas- 
sage distance r m \ n between the galaxies in Fig. [TO] it emerges 
that retrograde particles originate from further out in the in- 
falling disc. Prograde particles have average values of D of about 
Dp 1 ^ « 6 kpc, retrograde ones of approximately D^ ti0 « 7 kpc. 
For every model, the average value of D is higher for the retro- 
than the prograde particles. This effect is more obvious when 
only considering the far particles which have apocentre distances 
of at least 60 kpc. 

The rise in D with increasing r m ; n , by about 1 kpc from r m j n as 
8 to r m ; n w 16 kpc can be explained as follows: when the galaxies 
pass at greater distances, the tidal interaction is less severe. The 
number of stripped particles decreases with increasing r m ; n . As it 
is easiest to unbind particles from the rim of a galactic disc, these 



outside particles will still contribute to the pro- and retrograde 
populations, while particles with lower initial distances remain 
bound to the infalling galaxy. 

The orbital eccentricities for both particle populations rise 
slightly with increasing r m ; n . Increasing velocities have a 
stronger effect on the average orbit eccentricities for prograde 
particles, going from e pro = 0.74 to e pro = 0.97. The corre- 
sponding value for the retrograde population ranges only from 
<?retro = 0-79 to e retro = 0.90. As can be expected, far particles in 
all models show higher eccentricities on average. 

5.2. Equal-mass mergers 

For merging galaxies, both a pro- and a retrograde orientation of 
the infalling disc's direction of spin compared to the orbit lead 
to significant numbers of particles orbiting the merger remnant. 

In the case of merging galaxies, the outcoming distributions 
of pro- and retrograde particles seem not to follow a distinctive 
trend, the parameters do not depend in a simple way on the initial 
velocity or first pericentre between the two galaxies. The pro- 
grade orientation tends to give a higher absolute number of pro- 
and retrograde particles. N a \\ ranges from approximately 5000 
to 43000 for prograde infalling galaxies, it is between 6000 and 
27000 particles if the infalling galaxy spins retrogradely. For r m ; n 
above about 4 kpc, N a \\ is higher for a prograde than for a ret- 
rograde infalling galaxy, below it is smaller on average. These 
numbers drop when considering only far particles. The drop-off 
is most severe for models with low initial velocities (0.5 times 
the parabolic velocity), low r m i n (j 4 kpc) and a retrograde in- 
falling galaxy. Essentially no far particles are found for these. 
For faster and more distant encounters, several 10 3 to 10 4 par- 
ticles with high apocentres are counted. In principle, the maxi- 
mum distance of a satellite galaxy born as a TDG associated with 
a merger can be used to derive a lower boundary for the relative 
velocity between the precursor infalling galaxy and the central 
(target) galaxy. Therefore it can be concluded that mergers of 
galaxies with low relative velocities do not produce a sufficient 
number of far particles to resemble the MW satellite system. 

If the infalling galaxy's spin angular momentum is oriented 
retrogradely to the orbital one, f re tro varies between 0.79 to 1.00. 
When considering the far particles only, this fraction even drops 
to 0.38 in one model (2.5degl00vel), with an even spread in 
^retio- The retrograde fraction is plotted in Fig. [TT] against the 
first perigalactic passage distance of the galaxies, r m ; n . For a ret- 
rograde orientation, plotted in the right panels of Fig. [TT] frac- 
tions of retrograde particles vary between 0.01 and 0.37. Even 
^retro = 0-^0 is reached by one model when taking only the far 
particles into account. However, there is a large gap in the distri- 
bution, the values for the other models spread between 0.00 and 
0.25. 

The initial distances D within the infalling galaxy of where 
the tidal material stem from do not distinguish the pro- and ret- 
rograde populations as well as in case of fly-bys of galaxies. 
The average values differ in general, but the distributions over- 
lap strongly. Furthermore, there is no pronounced trend: for a 
retrograde spinning infalling galaxy, the average initial distance 
for prograde particles tends to be higher (at about D^ Q = 7 kpc) 

than that of retrograde ones (at about D^ tt0 = 6.5 kpc). This 
holds true up to a first perigalactic passage distance r m j n of about 
10 kpc. Thereafter, the three models with higher r m ; n show the 
opposite trend, with D^ Q of about 5 kpc. These differences are 
even smaller for far particles and above r m \ a = 10 kpc the far pro- 
grade particles even show slightly higher initial distances. For 
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Fig. 11. Equal-mass-mergers: the fraction of particles on retrograde orbits with apocentres above 30 kpc (P re ti , upper panels), and 
above 60 kpc (P l r f (m , lower panels) after the merger. The panels on the left show models with retrograde, those on the right prograde 
spin orientation of the infalling galaxy. See Tables lB~3"l to lB.6l for the values. The symbols reflect the initial velocities, with '050vel' 
standing for 0.5, 'lOOvel' for 1.0 and '150veF for 1.5 times the parabolic velocity, calculated assuming the initial galaxies are point 
masses. Mergers where the infalling disc has a retrograde angular momentum (left panels) lead to the highest fractions of retrograde 
material in the models computed. There is no simple dependence on the initial relative velocity of the galaxies, but for the subset of 
far particles, faster interactions tend to lead to higher values in P ret ro- Mergers where the infalling disc has a prograde spin angular 
momentum (right panels) lead to low fractions of retrograde material. Similar to the fly-by interactions, faster initial velocities lead 
to higher fractions of retrograde material. But in contrast to the fly-by cases, P retro does not drop when considering only the subset 
of far particles, because of the more similar radial distribution of particles on pro- and retrograde orbits. See the on-line edition of 
the article for a colour version of this figure. 



prograde infalling galaxies there is a slight trend to lower D for 
prograde than for retrograde particles. 



The average orbit eccentricities of the two populations dif- 
fer in general. For models in which the infalling galaxy has a 
retrograde spin, prograde particles show a higher average orbit 
eccentricity in most cases. This is the other way around in mod- 
els with prograde spinning infalling galaxies, there retrograde 
particles have higher eccentricities. These trends hold true when 
only considering far particles. The orbit eccentricities for these 
subpopulations increase, however, and the differences between 
pro- and retrograde particles are reduced. 



5.3. 4-to-1 mass fly-bys 

The results of the ten fly-by models with a mass ratio of 4: 1 for 
target:infalling galaxy are compiled in Tables |B~71 and lB~8l The 
spread in initial parameters here is higher, the initial velocity 
varies from 1.75 to 2.5 times the parabolic velocity, the peri- 
galactic distance of the two passing galaxies, r^, now ranges 
from 3.7 to 16.9 kpc. 

As in the equal-mass-situation, most models show a signifi- 
cant number of both pro- and retrograde particles, their sum be- 
ing in the range between 1300 and 14000. The fraction of ret- 
rograde particles P^L covers a wide range, from less than 0.01 
in model 8degl75vel up to 0.94 in model 4deg250vel. The num- 
bers of far particles fall in the range between 600 and 6000. As in 
the equal-mass models, P f r f tm is usually much lower, in the range 
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of several per cent. But for the fastest initial velocities, again, 
^retro does not change as much. Retrograde material again forms 
first and the transition of the two phases, producing at first retro- 
and then prograde particles, also takes place in interactions of 
non-equal-mass galaxies. It can therefore be concluded that the 
two-phase origin is a general behaviour of fly-by interactions, 
and not a special effect of a certain mass-ratio of the galaxies. 

One important difference to the equal-mass case is that P leao 
now seems to change dramatically with r m j n . For two models 
of the same initial velocity the retrograde-fraction can differ by 
a factor of 15 or more. Still, looking at the maximum values 
for Pretio, the previously found trend of an increasing fraction 
of retrograde material with increasing velocity is reproduced. 
Interestingly, these maximum values are all found for r m ; n in the 
short range of 6.9 to 9.7 kpc, r m ; n for the other models here are 
at least a factor 1 .5 off. This has to be contrasted to the equal- 
mass case where the relative differences in r^n are only 1.2 and 
the difference between maximum and minimum in r m \ n does not 
exceed a factor of two. 

It can be speculated that, if more calculations with interme- 
diate r m ; n would be included, the P retro -values would again lie on 
curves for each initial velocity that are not crossing each other, 
similar to Fig. [9] A reconstruction similar to the one described in 
Sect. lS.ll will be more difficult in this case because the fraction of 
retrograde material is not only determined by the initial velocity 
but also by the minimum distance to which the passing galaxy 
approached the target. Nevertheless, as probably every point in 
the f letro -r m i n -plane can be associated to a given initial velocity, 
such a reconstruction is deemed to be in principle possible in the 
case of a 4-to-l mass ratio. 

The average orbital eccentricities of the pro- and retrograde 
populations are similar to the equal-mass models. This includes 
the tendency that the retrograde population has a slightly lower 
orbital eccentricity in principle. The same is true for the site of 
origin in the precursor disc: again, retrograde particles originally 
had a larger radial distance D to the infalling disc's centre. The 
difference now on average is 0.75 kpc for all (af apo > 30 kpc) 
and 0.5 kpc for the far (c/ apo > 60 kpc) particles. For model 
2deg200vel no significant differences are found, while in others 
the difference can exceed 1.5 kpc. 

5.4. 4-to-1 mass mergers 

The tidal particle numbers in the merger models are again higher 
than in the fly-by case of this mass-ratio. The sum of pro- and 
retrograde particles N a \\ reaches up to 60000 particles and more 
in some models. The number of far particles reaches up to 30000 
for both, pro- and retrograde orientations of the infalling galaxy. 
While the numbers are somewhat higher than in the equal-mass 
case, one has to keep in mind that one particle now has only half 
the mass. Nevertheless, the total amount of material on pro- and 
retrograde orbits is comparable. 

In all 2 x 6 merger models, particles on pro- and on ret- 
rograde orbits are found. In contrast to the equal-mass merger 
models, the fraction of retrograde material, P retI0 , now is more 
evenly distributed. While in the equal-mass case the infall of a 
disc with retrograde spin leads to high values of P retro , in the 
4-to-l -mass case Pf^ tIO is in the intermediate regime of 0.23 to 
0.63. Considering only particles with apocentres above 60 kpc, 
^retro cnan g es to 0.03 and 0.06 for the two models with the 
biggest r m j n (they have the smallest Pj^-values, too: 0.23 and 
0.30 respectively). The fraction is not much different in the other 
four models. A prograde spin of the infalling galaxy results in 



lower values of P letm , but not as low as in the equal-mass case: 
PfX- ranges from 0.01 to 0.44, P f T f ao even reaches 0.66 in model 
2.5deg050vel. As far as the six models each tell, P letro drops with 
increasing r m ; n . The models with lower velocities have higher 
^retro, too, but are at the same time the ones with the lowest val- 
ues of r min . 

The average orbital eccentricities e are generally on the 
higher end of the distribution found for equal-mass mergers. 
They vary between e = 0.76 and e = 0.90, with the retrograde 
ones being higher on average. This holds true for both spin ori- 
entations of the infalling galaxy but is more pronounced for a 
prograde infalling galaxy. Considering the far particles only, the 
average orbital eccentricities are again higher. 

The initial distance of the tidal material to the centre of 
the infalling galaxy, D al1 , does not show a conclusive behaviour 
when looking at all particles. But, considering only the far parti- 
cles, D far is higher for retro- than for prograde particles (£>^ ro > 
Dp") in most models (the ones with faster initial velocities). This 
is the same trend already observed in all fly-by interactions, but 
interestingly not in the equal-mass mergers where, if at all, an 
opposite trend was evident. 

6. Discussion 

6.1. Observational Predictions 

To compare the results with observations, the ratio between pro- 
and retrograde particles can be expressed with respect to the 
lower-numbered population. This counter-rotating ratio /? COU ntrot 
is defined as the fraction of particles orbiting in an opposite 
sense compared to the majority of pro- or retrograde particles. 
By this definition any knowledge about the orbital direction of 
the galaxy-interaction is removed from the data derived from the 
models, just as is the case for the observable situation around 
galaxies. This implies the following definition: 

( jyretro . jypro jf jyretro < jypro^ 
^countrot — \ ^ypro . ^yietro j£ ^yretro ^ ^ypro 

6.1.1. Numbers, masses and radial distribution 

A range in /? C ountrot of less than 1 : 100 up to about 1 : 1 can 
be easily reproduced with fly-bys and mergers of both mass- 
ratios investigated. As all stellar particles in a calculation have 
the same mass, these number-ratios also reflect the mass-ratios. 
The maximum mass available can be calculated from the total 
number of orbiting particles. The mass of all particles orbiting 
the target galaxy with an apocenter larger than 30 kpc is up to 
6.8 x 10 8 Mq, 8.6 percent of the infalling galaxy mass in the 
equal-mass merger models. It is up to 9.2 x 10 7 Mq (1.2 per- 
cent) for equal-mass fly-bys, up to 5.6 x 10 8 Mq (14 percent) 
and 1.1 x 10 8 Mq (2.8 percent) for 4-to-l mass ratio mergers 
and fly-bys, respectively. As only particles that approach an apo- 
and pericentre are counted, these are only lower limits to the to- 
tal available mass. On average, a higher fraction of the infalling 
galaxy 's particles is stripped when it is less massive than the tar- 
get, favouring interactions of non-equalmass galaxies to disrupt 
the infalling galaxy more effectively. 

The two-phase origin of pro- and retrograde material in fly- 
by interactions results in a characteristic radial distribution of 
the orbits of tidal material. Particles on retrograde orbits are dis- 
tributed up to a maximum apocenter, while particles on prograde 
orbits can spread out much further. Thus, there is a maximum 
distance from the target galaxy up to where both co- and counter 
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orbiting material exists. Further out, only one orbital direction, 
the prograde one, can be found. This might be used as a hint on 
the original orbital direction when reconstructing past galaxy- 
galaxy interactions. 

6.1.2. Metallicity 

Even though the models do not include chemical evolution, some 
predictions on the metallicity distributions for pro- and retro- 
grade tidal material can still be deduced. It is a common fea- 
ture of disc galaxies to show a radial metallicity gradient, with 
metallicity decre asing with increasing rad ius D (see for exam- 
ple the review bv lHenrv & Worthevlll999l) . The infalling galaxy 
model has similar p arameters compared to those adopted by 
iMagrini et aT] d2007l) for M33. Therefore, their computed [Fe/H] 
gradient at more than 8 Gyr before today of -0.08 dex kpc 1 is 
taken for the model. Note that the following discussion strongly 
depends on the adopted metallicity gradient for the infalling disc 
galaxy. 

The different mean initial distances of the ultimately pro- 
and retrograde populations hint at a metallicity difference be- 
tween those particle groups. In the fly-by models, the mean ini- 
tial distance for prograde material is on average 0.9 kpc closer 
to the infalling galaxy's centre than that of retrograde material. 
This results in an average intrinsic [Fe/H] difference of 0.07 dex. 
In all fly-by models analysed, retrograde material has a lower 
metallicity compared to prograde material. This is irrespective 
of the mass ratio of the two galaxies. The intrinsic metallicity 
thus might be a useful tool to reconstruct the encounter-histories 
of stellar streams and maybe TDG systems around galaxies. 

The largest difference of average initial distances is 2.0 kpc 
in case of model 7deg200vel and only considering the 'far' par- 
ticles. This predicts a mean metallicity of prograde material of 
about 0.16 dex higher than that of retrograde material, given 
the adopted metallicity gradient in the infalling disc. The values 
strongly depend on the assumed metallicity gradient, but gener- 
ally lies b elow the uncertain ties of most metallicity determina- 
tions (e.g. lGrebel et al.f2003l) . It should become accessible once 
better statistics are available. However, the spread in initial dis- 
tance within one population is of the same order as the difference 
between the two. So, while a general trend can be expected, the 
metallicity can not be used to determine whether material or- 
bits pro- and retrograde. The results of merger models show that 
there the distribution in initial distance follows no clear trend. 

The scenario of a tidal origin for the MW streams and satel- 
lite galaxies of the MW predicts that in fly-by cases the intrin- 
sic metallicity of retrograde material is lower than that of pro- 
grade material. With better models including chemical evolu- 
tion and more precise observational data, it might become pos- 
sible to test this prediction: while the distributions of intrinsic 
metallicities of pro- and retrograde material overlap, retrograde 
material is expected to be slightly less metal rich on average.. 
But caution is advised: the accretion of unprocessed material 
by newly formed dwarf galax ies and the internal chemical self- 
evolution ( Rec chi et alj |2007) will most probably smear out this 
effect. It will thus be of importance to observe the initial (i.e. 
pre-encounter) metallicity as present in the oldest stars in stellar 
streams and dSphs. 

6.2. The MW satellites as TDGs 

The dSph satellite galaxies of the MW show properties that are 
difficult to understand in terms of them being cosmological sub- 



structures dMcGaugh & Wolfl2010h. Instead thev are more read- 
ily un derstood as ancient TDGs dMetz et al.ll2007t iRroupa et alj 
2010). Compared with the predictions of the presented models, 
it is found that the MW satellite system is also in agreement with 
having formed out of tidal material and a possible scenario for 
their origin can be distilled from the data. 

Is is sometimes observed that star-forming-regions in 
tidal tails are distr ibuted in a beads-on-a-string morphology 
(Mull anetalJl2011h . that is, star formation takes place at spots 
about equidistantly distributed along the tidal tails. This supports 
an approximately homogeneous distribution of the formation 
sites of TDGs along the t ails. Further hints for su ch a distribu- 
tion come from figure 2 in B ournaud & Dud (12006b . which plots 
the number of objects formed as a function of their formation 
radius, normalized by the total radial extend of the tidal tail. For 
their initial distribution (200 Myr after the interaction), the radial 
distribution is approximately constant, the bin-heights vary only 
by a factor of two. Therefore, in the following the most simple 
assumption on the formation efficiency of TDGs is adopted, their 
formed number being proportional to the number of particles on 
pro- and retrograde orbits. 

6.2.1. Orbital Poles 

The distribution of orbital poles of particles in models 
5deg200vel and 7.5degl00vel as presented in Figure [4] shows 
two distinct populations, offset by 180 degrees. This agrees with 
the positions of the MW satellite orbital poles, which show one 
concentration, from which Sculptor is offset by about 180 de- 
grees. In addition to that, it is especially interesting to note the 
shape of the distributions: they are elongated in diagonal direc- 
tions in both models. This is suggestive for the distribution of 
the orbital poles o f the co-orbiting MW satellites, as found by 
Metz et all d2008l) . Their figure 1 is reproduced in Fig. [T2l in 
the same projection as Fig.|4]here. Here the Zjyiw-axis is shifted 
to align with the positions of pro- and retrograde particles. The 
satellite galaxies Draco and Carina are oriented in a similar di- 
agonal direction compared to the other co-orbiting galaxies. 

The resemblance to the result for the merger model in Fig. [4] 
is most striking: not only are the most likely positions of the or- 
bital poles for the satellites reproduced by the model, but also the 
segments of the great circles (green lines), which represent the 
uncertainties of the orbital poles, align with the distribution of 
angular momenta. The counter-orbiting population is only made 
up of one galaxy, Sculptor, thus there is no structure that could 
be compared to the models. 

While the position of Sculptor's orbital pole is in agreement 
with the results of the tidal formation scenario assumed in this 
paper, further data on the proper motions of the MW satellite 
galaxies will allow the assessment of the validity of this sce- 
nario. Finding more orbital poles close to the one of Sculptor 
will allow to compare the angular-momenta substructures to the 
observed situation in more detail. This will test the validity of the 
tidal origin scenario and, if the MW satellite galaxies are TDGs, 
constrain the properties of the early interaction. 

6.2.2. Number- and mass-ratios 

Comparing the mass ratios and number ratios to the eight satel- 
lite galaxies with measured proper motions, the number-ratio 
of counter- to co-orbiting satellites in the known MW dwarfs 
is found to be ^countrot = 1:6 within the uncertainties, but 
^countrot = 1:4, excluding the Draco and Carina dwarf that show 
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Fig. 12. Possible orbital poles (dotted lines) in galactocentric galactic coordinates for the 24 known MW satellite galaxies as defined 
by their angular positions. For eight MW satellite galaxies proper motions are available, the resulting orbital poles are marked by 
the symbols. The green lines give the uncertainties in direction of the orbital poles. The red triangle marks the mean direction of the 
orbital poles of the six satellites close to Zmw = 180°, excluding Sagittarius and Sculptor. The red circle marks the spherical standard 
deviation from this average direction, the dashed circles are the regions 15° and 30° degrees fr om the normal of the DoS fitted to 
the 1 1 classical MW satellites. The data for the orbital poles are taken from Met z et aU(l2008l) . The possible orbital poles (dotted 
lines) are densest at /mw ~ 0° and 180°, where they are expected to lie in the case of a rotationally supported DoS. The known 
satellite orbital poles fall close to these positions with the exception of Sagittarius. A comparison of this plot with the results of 
particle orbital poles from the numerical calculations in Fig.|4]shows that the orbital poles of the MW satellite galaxies have a similar 
distribution. There are two distinct populations with an offset of about 180°. The positions of the Draco and Carina orbital poles 
even hint at a diagonal distribution of orbital poles (elongated from upper left to lower right), similar to the calculated distributions. 
This shows that the tidal origin model is in agreement with the observations of the MW satellite system. See the on-line edition of 
the article for a colour version of this figure. 



large orbital pole uncertainties, might be possible, too. These 
numbers lie well between the extremes of the models discussed 
in Sect. [4] 

The corresponding s tellar-mas s -rati os of the dwarf satellites 
are, using values from I Woo eta!] d2008l) : R™Zmt = 1 : 1700 
Considering the mass distributions, the model results (Sec. 16.1. H i 
have ratios of 1:100 to 1:1. But the observed ratio is strongly 
dominated by the LMC and SMC, which could be the infalling 
remnant of an initially passing galaxy (see Sect. 16.2.31 below). 
In that case, a ratio of about 1:17 excluding LMC and SMC is 
well reproduced by many of the far subpopulations, especially in 
the fly-by models with in the range between 0.05 and 0. 10 
and 0.95 to 0.99, but also for a number of merger models. 

Concerning the comparison of models with the MW, it is 
too early for a concluding statement because there are many 
dwarfs for which no proper motion measurements are avail- 
able. It is well possible that more dwarf galaxies on counter- 
rotating orbits will be found. Precise pro per-motion data will 
be available through the GAIA mission dPerrvman et al.l l2001: 
iLindegren et alj|2008l) , another method to determine proper mo- 
tions was presented by|Kaplinghat & StrigarJ d2008l) . 

Looking at the mass available in tidal material, the satel- 
lite dwarf galaxies of the MW have stellar masses in the range 
between a few 10 4 and 10 8 M Q , excepting LMC and SMC 
(Strig ariet al.ll2008l) . The maximum masses of pro- and retro- 



grade material in the calculations are sufficient to form several 
of such satellite galaxies, their number depending on the TDG 
formation rate in the tidal material. Judging from the two mod- 
els analysed in detail, the formation of retrograde tidal dwarf 
satellites of high mass seems to be unlikely. But in both calcula- 
tions there is plenty of mass on these orbits available to make up 
several Sculptor-type dwarfs. 

6.2.3. Possible reconstruction of the early MW encounter 

As an illustrative example, a possible early encounter the MW 
suffered shall now be reconstructed, under the assumption that 
the MW satellites formed from tidal material expelled during 
that interaction. As the results of merger models give no conclu- 
sive trends, the discussion in this section concentrates on fly-by 
models of equal mass for simplicity. 

The counter-rotating ratio /? C ountrot of the MW satellites is 
high, at least 1 : 4 considering satellite numbers only. Therefore 
in the models one has to look for P le t r o-values above 0.8 or be- 
low 0.2. As it is known that the satellite galaxies have apocentre 
distances above 60 kpc, only the results of far particles (lower 
panel in Fig. |9]l are considered. Several of the models might 
fit approximately and none can be excluded right away. Either 
the fly-by was fast (models with 2.4 times the parabolic veloc- 
ity), then fjetio is high, so what is seen today as counter-orbiting 
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satellite Sculptor was a prograde orbit compared to the relative 
motion of the interacting galaxies. Or P ret i is low, then all other 
models would fit and the counter-orbiting satellite comes from 
retrograde material. One way to discriminate between high or 
low Pretro would be the initial metallicity. 

The recent discovery of the most metal- poor star currently 
known in a MW satellite galaxy in Sculptor dFrebel et alj [2010) 
is in agreement with Sculptor having a lower average initial 
metallicity than other satellites. This might hint at an origin far- 
ther away from the centre of an infalling galaxy, which on aver- 
age is the case for retrograde material. Sculptor is, as a counter- 
orbiting satellite, in the minority, so P reao would thus have to be 
small. 

Reconstructing the initial encounter with these assumptions 
leads to the following scenario: about 10 Gyr ago the MW had 
an off-centre fly-by encounter with a similar disc galaxy which 
produced tidal dwarf galaxies. The infalling galaxy was on a 
parabolic to mildly hyperbolic orbit (not too fast as the models 
with 2.4 x Vpa ra b are excluded) with an orbital pole in the same 
direction as that of the later formed co-orbiting part of the DoS. 
But there remains a question: where did the infalling galaxy go'fl 

Maybe the initial interaction was not of similar-sized galax- 
ies but the infalling one was smaller and can indeed still be seen 
today: the Large Magellanic Cloud. Maybe it is now on its sec- 
ond approach to the MW. The orbital-pole-direction of the in- 
teraction must have been close to the normal vector of the Disc 
of Satellites, whic h points to /d s = 157.3° and Z?doS = -12.7° 
dMetz et al.ll2007l) . These values are similar to the orbital pole 
of the LMC, wh ich is / LMC = 172.0° and b hM c = -10.6° 
dMetz et all 120081) . The LMC spi n direction, how ever, can not 
be used to constrain this scenario. Weinberg (2000) predicts that 
the LMC disc precesses and nutates. T his was observationally 
supported by Ivan der Marel et al.l d2002l) . who found that the in- 
clination of the LMC disc currently changes at a rate of about 
100 degrees per Gyr. It is thus not possible to deduce the LMC's 
orientation several Gyr ago. The orbital direction of the LMC 
agrees with Sculptor being a retrograde satellite, as it is on a 
counter-rotating orbit, leading to a low P le tro- This is produced by 
most 1 : 1 and 4: 1 models, especially at larger distances from the 
target galaxy. If an extrapolation of the results to more extreme 
mass-ratios is possible, an assumption which is supported by the 
finding that less massive infalling galaxies are more effectively 
disrupted, the results are thus compatible with this scenaricQ. A 
low P le tro then hints at a not too fast fly-by (the highest-velocity 
models can be excluded, see Fig. [9]), in better agreement with 
an orbit that does not allow the LMC to escape. Furthermore, 
most models lie in the regime of lower P re tro values, therefore 
no fine-tuning of the initial parameters is needed to arrive at 
situations similar to today 's MW satellite system. Additionally, 
fly-by models with lower initial velocity lead to higher parti- 
cle numbers stripped from the infalling galaxy, thus more ma- 
terial is available to form several TDGs. The radial distribution 
of the observed satellites supports this, too: Sculptor is not too 
far from the MW 80 kpc) and has one of the lowest radial 
velocities, suggesting a smaller apocentre distance than for the 
other satellites. This is the natural outcome for retrograde parti- 
cles due to the phase-one-origin in the fly-by models, as can be 



7 Assuming a relative velocity of about 250 km/s, after 10 Gyr it will 
have a distance of 2.5 Mpc, at the fringes of the Local Group. 

8 However, a differing accretion history of the MW and the LMC 
progenitor might have made the MW grow more rapidly, thus the initial 
masses might have been more similar. 



seen in Fig.|7]for model 5deg200vel. That such models naturally 
form a disc of tidal debris is demonstrated in Fig. [5] 

This work has demonstrated above all else that the orbits of 
stellar streams and the dSph satellites will be a powerful tool for 
reproducing the early events that shaped the MW if they are in- 
terpreted to be of tidal origin. A more detailed modelling of this 
scenario is deemed to be of huge interest for the search of the 
origin of the MW satellite distribution. This would be in beau- 
tiful agreement with Lynden-Bell's (1 19761) original suggestion 
that the LMC could be the major surviving part of a Greater 
Magellanic Galaxy, parts of which have been torn off by tides 
forming dSph galaxies and streams. 

7. Conclusions 

In an attempt to approac h a solution to the Fritz-Zwicky para- 
dox dRroupaet al. 2010) and the peculiar spacial arrangement 
and internal properties of the Milky Way satellite galaxies, the 
generic properties of tidally expelled material in galaxy-galaxy 
interactions at an early cosmological epoch are here studied nu- 
merically. The galaxy models are set up by scaling the MW prop- 
erties down to those of an early, bulge-less galaxy at z ~ 2, about 
10 Gyr ago. 

Computations of interactions between two similar disc 
galaxies as well as those with mass-ratios of 4-to-l for infalling 
to target galaxy have been performed. It has be shown for the 
first time that it is possible to form particles in tidal streams with 
both pro- and retrograde orbits in one single galactic encounter. 
This is true for mergers and, maybe more surprisingly, for fly- 
bys, of nearly perpendicularly oriented disc galaxies. The typical 
origins for pro- and retrograde material are sketched in Fig. [13] 
In the merger case, the close passing of the infalling galaxy near 
the centre of mass of the target results in particles passing the 
centre on both sides. They therefore end up on opposing orbits 
around the merged galaxies. In addition, material tidally slung 
away falls back at the merged galaxies in pro- and retrograde 
directions. In the case of a fly-by, the galaxies depart and the 
infalling one evolves a tidal tail, from which particles fall back 
towards the target galaxy. The direction of the orbital motion for 
particles falling in is changed from retrograde, with respect to 
the interaction geometry, to prograde, after the tidal tail sweeps- 
over the centre of the target galaxy. This marks two phases of the 
encounter, effectively forming two classes of particles on oppos- 
ing orbits in a fly-by encounter. 

TDGs will form from the tidal material, occupying the same 
phase-space volume because the effect of dynamical friction on 
them can be neglected. This scenario explains the otherwise un- 
likely distribution of satellite positions in a Disc of Satellites 
and orbital poles and naturally resolves the riddle of the counter- 
orbiting satellite Sculptor. The number of particles on pro- and 
retrograde orbits are in agreement with common satellite galaxy 
masses of some 10 7 Mq, up to 14 per cent of the infalling 
galaxy's particles finally are on orbits around the target galaxy 
with apocentre distances beyond 30 kpc. 

Moreover, it is possible to extract constraints on the early 
properties of such a disc of tidal material. While the situation 
in case of a galaxy merger is complicated and depends strongly 
on the actual parameters of the interaction (initial relative ve- 
locity, perigalactic passing distance), fly-by events show a clear 
trend: the fraction of retrograde particles varies strongly with 
the initial relative velocity between the two galaxies, increasing 
for higher velocities. Thus, the fraction of retrograde material 
observed constrains the relative velocity at which the encounter 
happened. Considering particles with an apocentre of 60 kpc or 
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Fig. 13. Sketch of the typical origins for pro- and retrograde ma- 
terial in merger (top) and fly-by (bottom) interactions. In both 
cases, the infalling galaxy starts on the right and follows the path 
sketched by the dashed line. The direction of spin of its disc is 
indicated by the arrows along the edge in the fly-by sketch (pro- 
grade with respect to the orbit), for mergers both orientations 
of spin were modelled. The target galaxy is seen edge on. In a 
merger, the collision of the two galaxies involves particles of the 
infalling galaxy to pass the centre of mass of the target on both 
sides, resulting in pro- and retrograde orbits. As a merger might 
involve several passings of the galaxies before the final merging, 
other sources of pro- and retrograde material, like tidal tails, also 
play a role. In a fly-by, after the galaxies pass, the infalling one 
developes tidal tails. The tidal material streaming towards the 
target galaxy first follows back the infalling's path, resulting in 
retrograde material appearing first (Phase 1). Later, as the galax- 
ies further depart from each other, the tidal tail is swept over the 
target's centre of mass. This changes the orbital direction of the 
particles falling towards the target galaxy, making their orbits 
prograde (Phase 2). 



more, the increase is not as strong as retrograde particles have 
lower apocentres. Prograde particles fall in later when the galax- 
ies are more separated again, giving them higher apocentre dis- 
tances. They are also found to be on more eccentric orbits ini- 
tially, but so me circularisation of th eir orbit can be expected (see 
for example Hashimoto et al. 20031). 

Of observational interest is the site of origin of pro- and 
retrograde particles in the infalling galaxy's disc in fly-by en- 
counters. Prograde ones originate from a more central position 
by about 1 kpc. Combined with a metallicity gradient in the in- 
falling disc this results in an average difference in metallicity of 
pro- and retrograde material, with the latter being less metal rich 
by about 0.1 dex in [Fe/H]. The effect, however, will have been 
strongly affected and probably been washed out l ater because of 
individual chemical evolution (iRecchi et alj|2007l) and accretion 
of unprocessed material from the inter-galactic medium onto 
condensed objects such as TDGs. The results of merger mod- 
els are not as conclusive concerning the metallicities of pro- and 
retrograde material. 



Within the scenario of a tidal origin of the MW satellite 
galaxies, some first constraints on the early, major encounter 
the MW experienced can be derived from the present results. A 
merger with a small initial relative velocity of the two galaxies 
can be excluded as the origin of the MW dwarf satellite galaxies 
because no particles with high enough apocentre distances are 
produced. In principle, the maximum distance of a MW satel- 
lite born as a TDG associated with a merger can be used to 
derive a lower boundary for the relative velocity between the 
merging galaxies. In a fly-by scenario, an early interaction with 
the LMC progenitor shows an interesting consistency: the ini- 
tial direction of motion of the infalling galaxy is known, making 
Sculptor retrograde. This results in a low fraction of retrograde 
material, in agreement with the results of the majority of models. 
Furthermore, this fits with Sculptor's comparably lower apocen- 
tre and predicts a lower initial metallicity, as discussed above0. 
Such a scenario is in beautiful agreement with Lynden-Bell's 
(11976b original suggestion that a Greater Magellanic Galaxy 
might be the parent of the then-known dSph satellite galaxies. 
It might be the most probable origin of the MW satellite distri- 
bution, but needs to be addressed in more detail in the future. 
The implications of a tidal origin of th e MW satellite galax ies 
for cosmological theory are expanded in Kroupa et al.l d2010t) . 

The models discussed in this contribution did not include gas 
dynamics, nor was the growth of the Milky Way disc after the in- 
teraction included. No direct formation of satellite galaxies was 
investigated. These points will be addressed and included in fu- 
ture works, but will mainly change the detailed parameters of 
pro- and retrograde particle populations, not the fact that both 
orbital directions are produced naturally in interactions of disc 
galaxies. 

Acknowledgements. M.P. acknowledges support through the Bonn-Cologne 
Graduate School of Physics and Astronomy and through DFG reseach grant 
KR 1635/18-1 in the frame of the DFG Priority Programme 1177, Witnesses 
of Cosmic History: Formation and evolution of galaxies, black holes, and their 
environment. We thank Manuel Metz and Ingo Thies for very useful discussions. 



References 

Babusiaux, C, Gomez, A., Hill, V., et al. 2010, A&A, 519, A77+ 

Ballero, S. K., Matteucci, F., Origlia, L., & Rich, R. M. 2007, A&A, 467, 123 

Binney, J. & Tremaine, S. 1987, Galactic dynamics (Princeton, NJ, Princeton 

University Press, 1987, 747 p.) 
Boily, C. M., Kroupa, P., & Penarrubia-Garrido, J. 2001, New Astronomy, 6, 27 
Bournaud, F, Bois, M., Emsellem, E., & Due, P. 2008a, Astronomische 

Nachrichten, 329, 1025 
Bournaud, F. & Due, P. -A. 2006, A&A, 456, 481 
Bournaud, F., Due, P. -A., & Emsellem, E. 2008b, MNRAS, 389, L8 
Diemand, J., Kuhlen, M., Madau, P., et al. 2008, Nature, 454, 735 
D'Onghia, E. & Lake, G. 2009, in IAU Symposium, Vol. 256, The Magellanic 

Group and the Seven Dwarfs, ed. J. T. van Loon & J. M. Oliveira, 473 
Fellhauer, M. & Kroupa, P. 2002, MNRAS, 330, 642 

Fellhauer, M., Kroupa, P., Baumgardt, H., et al. 2000, New Astronomy, 5, 305 

Frebel, A., Kirby, E. N., & Simon, J. D. 2010, Nature, 464, 72 

Gerhard, O. 2002, Space Science Reviews, 100, 129 

Gilmore, G, Wilkinson, M. I., Wyse, R. F. G., et al. 2007, ApJ, 663, 948 

Grebel, E. K., Gallagher, III, J. S., & Harbeck, D. 2003, AJ, 125, 1926 

Hartwick, F. D. A. 2000, AJ, 119, 2248 

Hashimoto, Y., Funato, Y., & Makino, J. 2003, ApJ, 582, 196 

Henry, R. B. C. & Worthey, G. 1999, PASP, 111, 919 

Hemquist, L. 1990, ApJ, 356, 359 

Hemquist, L. 1993, ApJS, 86, 389 

Hunsberger, S. D., Charlton, J. C, & Zaritsky, D. 1996, ApJ, 462, 50 
Kaplinghat, M. & Strigari, L. E. 2008, ApJ, 682, L93 
Kerr, F. J. & Lynden-Bell, D. 1986, MNRAS, 221, 1023 

9 Note that the LMC spin is likely to be precessing and nutating, 
such that its current direction can not be us ed to constrain this scenario 
(Weinberg 200CMvan der Marel et alj|2002l) . 



19 



Pawlowski, Kroupa & de Boer: Making Counter-Orbiting Tidal Debris 



Klessen, R. S. & Kroupa, P. 1998, ApJ, 498, 143 

Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, ApJ, 522, 82 
Kroupa, P. 1997, New Astronomy, 2, 139 

Kroupa, P., Famaey, B., de Boer, K. S., et al. 2010, A&A, 523, A32+ 
Kroupa, P., Theis, C, & Boily, C. M. 2005, A&A, 431, 517 
Kunkel, W. E. 1979, ApJ, 228, 718 
Li, Y.-S. & Helmi, A. 2008, MNRAS, 385, 1365 
Libeskind, N. I., Frenk, C. S., Cole, S., et al. 2005, MNRAS, 363, 146 
Libeskind, N. I., Frenk, C. S., Cole, S., Jenkins, A., & Helly, J. C. 2009, 
MNRAS, 399, 550 

Lindegren, L., Babusiaux, C, Bailer-Jones, C, et al. 2008, in IAU Symposium, 
Vol. 248, The Gaia mission: science, organization and present status, ed. W. J. 
Jin, I. Platais, & M. A. C. Perryman, 217 

Lopez-Corredoira, M., Cabrera-Lavers, A., Garzon, F., & Hammersley, P. L. 
2002, A&A, 394, 883 

Lynden-Bell, D. 1976, MNRAS, 174, 695 

Lynden-Bell, D. & Lynden-Bell, R. M. 1995, MNRAS, 275, 429 
Madejsky, R. & Bien, R. 1993, A&A, 280, 383 
Magrini, L., Corbelli, E., & Galli, D. 2007, A&A, 470, 843 
Majewski, S. R. 1994, ApJ, 431, L17 

Mao, S., Mo, H. J., & White, S. D. M. 1998, MNRAS, 297, L71 

McGaugh, S. S. & Wolf, J. 2010, ApJ, 722, 248 

Metz, M., Kroupa, P., & Jerjen, H. 2007, MNRAS, 374, 1 125 

Metz, M., Kroupa, P., & Jerjen, H. 2009, MNRAS 

Metz, M., Kroupa, P., & Libeskind, N. I. 2008, ApJ, 680, 287 

Metz, M., Kroupa, P., Theis, C, Hensler, G., & Jerjen, H. 2009, ApJ, 697, 269 

Minniti, D. & Zoccali, M. 2008, in IAU Symposium, Vol. 245, The Galactic 

bulge: a review, 323 
Mirabel, I. E, Dottori, H., & Lutz, D. 1992, A&A, 256, L19 
Moore, B., Ghigna, S., Govemato, E, et al. 1999, ApJ, 524, L19 
Mullan, B., Konstantopoulos, I. S., Kepley, A. A., et al. 2011, ArXiv e-prints 
Naab, T. & Ostriker, J. P. 2006, MNRAS, 366, 899 
Ojha, D. K. 2001, MNRAS, 322, 426 
Okazaki, T. & Taniguchi, Y. 2000, ApJ, 543, 149 
Palma, C, Majewski, S. R., & Johnston, K. V. 2002, ApJ, 564, 736 
Perryman, M. A. C, de Boer, K. S., Gilmore, G., et al. 2001, A&A, 369, 339 
Recchi, S., Theis, C, Kroupa, P., & Hensler, G. 2007, A&A, 470, L5 
Sales, L. V, Navarro, J. E, Abadi, M. G.. & Steinmetz, M. 2007, MNRAS, 379, 

1464 

Sawa, T. & Fujimoto, M. 2005, PASJ, 57, 429 

Spergel, D. N, Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175 

Stoehr, E, White, S. D. M., Tormen, G., & Springel, V. 2002, MNRAS, 335, L84 

Strigari, L. E., Bullock, J. S., Kaplinghat, M., et al. 2007, ApJ, 669, 676 

Strigari, L. E., Bullock, J. S., Kaplinghat, M., et al. 2008, Nature, 454, 1096 

Tollerud, E. J., Bullock, J. S., Strigari, L. E., & Willman, B. 2008, ApJ, 688, 277 

Tully, R. B. & Fisher, J. R. 1977, A&A, 54, 661 

van der Kruit, P. C. 1987, A&A, 173, 59 

van der Marel, R. P., Alves, D. R., Hardy, E., & Suntzeff, N. B. 2002, AJ, 124, 
2639 

Vogt, N. P., Phillips, A. C, Faber, S. M., et al. 1997, ApJ, 479, L121 
Walter, E, Martin, C. L., & Ott, J. 2006, AJ, 132, 2289 

Weilbacher, P. M., Due, P. -A., & Fritze-v. Alvensleben, U. 2003, A&A, 397, 545 
Weinberg, M. D. 2000, ApJ, 532, 922 

Wetzstein, M., Naab, T, & Burkert, A. 2007, MNRAS, 375, 805 
White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341 
Woo, J., Courteau, S., & Dekel, A. 2008, MNRAS, 390, 1453 
Zhao, H. 1998, ApJ, 500, L149 
Zwicky, F. 1937, ApJ, 86, 217 

Zwicky, F. 1956, Ergebnisse der exakten Naturwissenschaften, 29, 344 



Appendix A: Scaling of the Galaxy Models 

A.1. Assumptions 

The model is based on several assumptions used to determine 
the main structural properties of the disc galaxy model by scal- 
ing current MW parameters to a smaller disc mass. A constant 
and time-independent proportionality between luminosity and 
mass is assumed, to allow to treat surface brightness and surface 
density similarly. Th is allows to apply the Tully-Fisher relation 
(Tully & Fisher l 19771) to scale the circular velocities. Thus, also 
the validity of the Tully-Fisher relation independent of age or 
redshift and size of the disc galaxy is assumed, whic h seems to 
be acceptable, at least to a redshift z ~ 1, following Vog t et all 



Table A.l. Assumed MW parameters. 



Symbol 


Value 


Description 


M ,oaa y 


5.6 x 10 10 M o 


mass of disc 


•-.today 

d , 


1.3 to 1.8 x 10 10 M© 


mass of bulge 


3.5 kpc 


disc scale length 


(-.today 

R Q 


8.5 kpc 


Solar radius 


today 
*max 


220 km s -1 


circular velocity near Sun 




50 M pc~ 2 


disc surface density near Sun 



Assumed parameters of today's MW disc, applied as starting points for 
the down-scaling of the disc to its younger version, 10 Gyr ago. 



(119971) . iMaoet al. (1998) have found no strong indications for 
a change in the relation, too. The gaseous galaxy-component is 
ignored, assuming a progenitor galaxy described by a pure expo- 
nential disc embedded in a spherical Hernquist dark-matter halo, 
which is used to model the empirically found fiat rotation curves 
in disc galaxies and constitutes, in a non-Newtonian framework, 
the zeroth-order approximation to non-Newtonian dynamics. 

A.2. Determination of parameters 

To determine the parameters of the model, the stellar mass of the 
bulge component of the MW as seen today is taken as a start- 
ing p oint. Values range from 1.6 x 10 10 M0 dMinniti & Zoccahl 
2008 ) to 1.8 x 10 10 M r) (calculated using values from 
iBinnev & Tremaind[T987l) or might be as low as 1 .3 x 10 10 M Q 



(lGerhardll2002ir Based on the presupposition that a merger of 
two identical galaxies created the bulk of the bulge component, 
the mass of the model galaxy has to be about half of today's 
bulge mass. Some mass loss might occur due to matter trans- 
ported outwards along tidal arms. This leads to a stellar mass of 
Md = 8.0 x 10 y Mq for the initial small disc galaxies in the case 
of the equal mass encounters. For the 4-to-l mass ratio interac- 
tions, disc masses of 4.0 x 10 9 Mq and M d = 1.6 x 10 10 Mq 
are adopted. After formation of the bulge the early Milky Way 
would still have accreted matter in a disc, forming today's disc 
component within about 10 10 years when accreting ~ 5 Mq yr~' . 
The dark matter halo component is chosen to be ten times more 
massive than the early disc. 

Furthermore, the following common values are adopted as 
those of today's MW and applied as starting points for the down- 
scaling: the disc scale length is chosen to be /? t d ° ddy = 3.5 kpc. 
Observations reveal similar values like 3.3 kpc from star c ounts 
from 2MASS performed by Lopez-Corredoira et al. (2002), but 
lower values, like 2.1 kpc (Gerhard 2002), have been reported, 
too. The larger value is chosen here because the thick disc seems 
to have a larger scale length (3.7 kpc according to Ojha 2001) 
and, as it is composed of older stars, should trace the original 
mass distribution better. It is a ssumed that the mass of the MW 
disc is M d oday = 5.6 x 10 10 M g . lGerhardld2002l) gives a value for 
the mass enclosed within 8 kpc of 5.5 x 10 10 Mq including the 
bulge . Also, the International A stronomical Union (IAU) stan- 
dard (Ke rr & Lynden-Bell 1 1986b for the circular velocity near 



the Sun v, 
the Sun R 



today 
max 
today 



= 220 km s and the Galactocentric distance of 



- 8.5 kpc is adopted. The adopted parameters are 
compiled in Table |A.f1 
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A.3. Estimating the size of the young MW 

The most important parameters to determine for an exponential 
disc are the scale length and the rotation curve. Furthermore the 
disc should be stable, even though small galaxies are often ob- 
served to be irregular (e.g. the LMC). 

The circular velocity v mt {R) at a radius R around a mass dis- 
tribution characterised by an enclosed mass of M(< R) is given 
by 



G ■ M{< R) 



R 



(A.l) 



where G is the gravitational constant. Using a two-component 
galaxy, the rotation curve is a result of the interplay between 
the disc- and the halo-mass-distributio n, the respectiv e circular 
velocities being added in quadrature dHernquisC l 99 3 ) . 

As a first estimate for the disc scale length, the central 
surface brightness is demanded to be cons tant for all galaxies 
and inde pendent of the Hubble par ameter ( Ivan der Kruitl [l987. 
see also iBinnev & Tremainel[r987f) . With the previously men- 
tioned assumption of a proportionality between luminosity and 
mass, this translates int o a constant central surface density "Eq. 
Naab & Ostriker (2006) applied the same assumption of a con- 
stant central surface density when studying the evolution of the 
MW. For the central surface brightness the following equation 
holds in the case of a thin exponential disc with a total mass Md 
and a scale length R d , 



2n = 



M d 



(2n-R 2 d )' 



(A.2) 



This leads to the scaling relation 



larger scale length of 1 .6 kpc is choosen for the model, because 
a bigger value increases the stability of the disc. This is because 
more halo mass (in relation to disc mass), with particles of ran- 
dom motions, is inside the disc's volume. The disc is truncated 
at 5 times this scale length, clo se to the frequently used value of 
4.5 times (Ivan der Kruitlll987l) . Note that the disc of a gas com- 
ponent today is understood to be more extended than the stellar 
disc. 

The rotation curve is demanded to be flat and the maxi- 
mum velocity is estimated by applying the Tully-Fisher relation. 
The maximum rotational velocity to adjust the density profiles 
to is thus scaled using the mass ratio of today's and the model 
galaxy's disc mass. This leads to: 

vj J = 115. ..135 km s _I , 

depending on the slope of the Tully-Fisher relation. An interme- 
diate value is adopted, aiming at a rotation velocity of v tot « 
125 km s" 1 . 

As the primary objective is to set up the disc component, 
whose mass-distribution is already fixed with the scale length, 
the rotation curve is adjusted with the dark matter halo. The halo 
properties are thus determined by demanding it to give that cir- 
cular velocity. 

For the disc, the rotational velocity v dmc at the cut-off radius 

t— cue 

can be estimated using Equation IA.1I inserting the total disc 
mass Md and the maximum disc radius 7?™ x . Having v tot and 

u d & cue 

y disc y haio j s eas jjy calculated making use of the fact that circular 

circ ' circ J ° 

velocities add in quadrature. This results in v hal ° = 105 km s . 

j — etfc ■ — i — i 

In a Hernquist dark-matter halo profile (Hernquist 1990), the 

density p n at a radial distance R is given by 



R 



today 



M 6 



ph = 



today ' 
d 



2n R(R + fl) 3 



Such that R d = 1.3 kpc. 

Another possibility to scale R d avoids inclusion of the MW 
disc mass but starts from the disc surface density S(/?Q ddy ) in the M{< R) = Mh 
Solar neighbourhood. Assuming an exponential disc, the central 
MW surface density Eo is calculated via 

Z(^ day ) = Z exp(-^ day /< day ) 



with the core radius a and the total halo mass Mh- The mass 
inside a radius R is 



(A.3) 



and then Eq. IA.2I is solved for the scale length, yielding R& = 
1.5 kpc when accepting a value of £(^|^ day ) = 50 Mq pc~ 2 as 
indicated from observations (see Naab & Ostriker 2006, and ref- 
erences therein). 

An estimate completely independent of today's MW prop- 
erties can be derived fro m typical properties of galaxies. 
Binn ev & Tremainel d 1987b stated that the typical disc scale 
length for galaxies scales with R& = 3/h kpc, where h 
is the dimensionless Hubble-parameter defined by Ho = 
100 h km s Mpc~' . Assuming a closed universe with currently 
Ho = 71 km s" 1 Mpc , mass-density parameter Qm = 0.27 
from Spergel et al. (2003) and vacuum density Q vac = 0.73, the 
time about 10 Gyr ago corresponds to a redshift of z ~ 2. The 
Hubble parameter at that redshift is given by 

H(z) = H [o vac + (1 - Q vac - Qm)(1 + zf + Q M (1 + z) 3 ] ' 5 , 

and was approximately H(z — 2) « 200 km s _1 Mpc 1 at that 
time, so h(z = 2) « 2.0. This again leads to R d — 1 .5 kpc which is 
in good agreement with the first two estimates. Finally, a slightly 



Again making use of Equation lA.il the mass within R = /?™ x , 
the maximum disc particle radius, is determined. With the use of 
Eq. lA.3l the halo core radius can be estimated. 

The model's halo needs to have 1 .9x 10 10 M0 within a radius 
of /?™ x = 8.0 kpc. In the case of a finite Hernquist profile with 
maximum radius 7?™" = 100 kpc this results in a core radius of 
a = 10.0 kpc. After virialisation, the rotation curve of the disc 
component has a maximum at * 122 km s . 

It should be emphasized that this 'dark-matter' halo has not 
been constructed for consistency with the CDM hypothesis, but 
merely to model the actual gravitational behaviour of the em- 
pirical flat rotation curves of the young galaxy. What the true 
explanation for this grav itational behaviour might be is the topic 
of other studies (see e.g. Kroupa et al. 20101) . 

For a quicker reference, the model parameters are compiled 
in Table Q] Interestingly, the resulting values are quite similar to 
those fou nd for M33, which is a bulge-less galaxy, too (see for 
example Mag rini et alj |2007). justifying the choice of parame- 
ters as physically reasonable. 

Furthermore, several models with a mass ratio of 4:1 be- 
tween the target and the infalling galaxy are calculated. Their 
parameters are determined in the same manner and can also be 
found in Table[l] while their results are discussed in Sect. [5] 
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A.4. Model setup 

The model gala xy is created using the code MaGalie 
dBoilv et al.ll2001l) . It is based on a method by iHernquisH d 1993b 
to generate multi-component galaxies. The code is highly effi- 
cient for generating multi-component disc galaxies with many 
particles as the computational time scales only linearly with the 
particle number. 

The disc components are set up with 5 x 10 5 particles, the 
halo components are made up of 10 6 particles. These numbers 
represent a compromise between resolution (demanding high 
particle numbers) and calculation time. They result in disc par- 
ticle masses of 1.6 x 10 3 M© for the equal-mass model or 
0.8 x 10 3 M© and 3.2x 10 3 M© for the infalling and target galaxy 
in the 4-to-l mass ratio models. 

As the integration code Superbox++ demands all particles in 
one galaxy to be of equal mass, the following standard work- 
around is used to allow the modelling of the stellar disc with 
a suitable number of particles without increasing the total par- 
ticle number to an unfeasible value. The idea is to split up the 
galaxy into a disc and a halo component and to treat them as 
separate galaxies whose centres of mass have the same posi- 
tion. That way, the mass per halo particle can be higher, result- 
ing in fewer halo particles and thus speeding up the computa- 
tions. Furthermore, the splitting in two models allows to tailor 
the grid sizes to the respective components' size and therefore 
gives a higher resolution to the important disc components. The 
total number of grids is doubled by this splitting, leading to an 
increase in computational time by about a factor two. But the 
reduced computation time through the smaller number of halo 
particles, scaling about linearly, outweights this. 

The number of grid cells per dimension is 64. The gridsizes 
(radii) are 10, 30, and 400 kpc for the disc component galaxy 
model and 40, 120 and 400 kpc for the dark matter halo com- 
ponent. This gives grid cells with side lengths of 0.33, 1.0 and 
13.33 kpc for the disc component and 1.33, 4.0 and 13.33 kpc 
for the halo component. 

The newly set up model still suffers from minor numerical 
set-up mismatches because velocities and positions are stored as 
discrete values. Consequently the galaxy has to virialise to arrive 
at dynamical equilibrium. For this the model, consisting of disc 
and halo component, is calculated in isolation with Superbox++ 
for 3000 steps with a time-step size of 0.5 Myr, corresponding 
to about 7 orbiting times at a radius of 4 kpc. The galaxy is in 
dynamical equilibrium thereafter. The total energy error during 
virialisation is less than 0.2 per cent, the Lagrange-Radii expand 
slightly as the sharp edge at the maximum particle distance cut- 
off is smoothed out. 

Appendix B: Tables of parameter scans 
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Table B.l. Equal-mass fly-by: particle numbers 



Model 


^ini [^parab] 


r- [kpc] 


Steps analysed 




jyfiir 


^retro 


^retro 


jwoU— zjUkpc 


D 60-z50kpc 

retro 


5degl80vel 


1.8 


7.44 


15 


4898 


805 


0.19 


0.01 


1641 


0.43 


6degl80vel 


1.8 


9.39 


15 


2874 


474 


0.13 


0.00 


281 


0.17 


7degl80vel 


1.8 


11.52 


15 


2011 


387 


0.15 


0.00 


322 


0.07 


8degl80vel 


1.8 


13.81 


15 


1476 


428 


0.09 


0.00 


362 


0.00 


5deg200vel 
6deg200vel 


2.0 
2.0 


8.35 
10.54 


15 
15 


5747 
4644 


2090 
1316 


0.40 
0.44 


0.08 
0.07 


1741 
1119 


0.02 
0.04 


7deg200vel 
8deg200vel 


2.0 
2.0 


12.91 
15.45 


15 
15 


2848 
1636 


902 
655 


0.43 
0.42 


0.05 
0.01 


663 
583 


0.02 
0.01 


5deg220vel 


2.2 


9.15 


15 


3233 


1372 


0.56 


0.13 


1292 


0.06 


6deg220vel 


2.2 


11.57 


15 


3836 


1195 


0.63 


0.09 


1046 


0.03 


7deg220vel 


2.2 


14.15 


15 


2177 


650 


0.66 


0.10 


616 


0.06 


8deg220vel 


2.2 


16.89 


15 


849 


364 


0.51 


0.09 


393 


0.04 


5deg240vel 


2.4 


9.98 


15 


576 


405 


0.73 


0.62 


694 


0.28 


6deg240vel 
7deg240vel 


2.4 
2.4 


12.56 
15.29 


15 
15 


2425 
1873 


1173 
693 


0.93 
0.98 


0.85 
0.95 


1171 

607 


0.55 
0.63 


8deg240vel 


2.4 


18.17 


15 


862 


479 


0.99 


0.99 


363 


0.80 



Vm: initial relative velocity of the two galaxies in units of the point-mass parabolic velocity. 
r min : pericentre distance between the two galaxy centres of density for the first passage. 

Steps analysed: number of snapshot files (in steps of 0.5 Gyr) analysed for determining pro- and retrograde motions. Particles had to appear 
within the circles of acceptance in 12 out of these snapshots. In those cases where the galaxies merged late, fewer than 12 snapshots have been 
analysed. In these cases, the particles had to appear within the circles in all snapshots. 
N: sum of pro- and retrograde particles. 

P rctro : fraction of retrograde particles with respect to the sum of pro- and retrograde ones. 

iV60-250kpc. sum of pro- and retrograde particles with galactic distances in the range of 60 to 250 kpc at the final snapshot. 
^re°trc 250kpC : fraction of retrograde particles in N 6 °- 2S0k ^. 

Table B.2. Equal-mass fly-by: properties of the two populations 



Model gr C ro gg e% m ZCtkpc] jgLgg^ gPqg g|o [kpc] 

5degl80vel 0.74 ±0.11 0.79 ± 0.12 0.75 ± 0.10 0.85 ± 0.05 5.67 ± 1.19 7.26 ± 1.16 5.97 ± 1.34 5.09 ± 0.39 

6degl80vel 0.78 ± 0.09 0.81 ±0.11 0.82 ± 0.07 - 5.71 ± 1.17 7.16 ± 1.25 5.83 ± 1.22 

7degl80vel 0.82 ± 0.08 0.82 ± 0.08 0.81 ± 0.06 5.91 ± 1.23 7.02 ± 1.36 5.95 ± 1.31 

8degl80vel 0.84 ± 0.06 0.84 ± 0.07 0.85 ± 0.04 - 6.01 ± 1.24 7.55 ± 1.21 5.98 ± 1.22 - 

5deg200vel 0.82 ± 0.08 0.80 ±0.11 0.80 ± 0.08 0.93 ± 0.03 5.89 ± 1.05 6.73 ± 0.68 6.07 ± 1.07 7.23 ± 0.71 

6deg200vel 0.83 ± 0.06 0.79 ±0.12 0.84 ± 0.04 0.90 ± 0.04 5.71 ± 1.05 6.65 ± 1.17 5.85 ±1.11 7.78 ± 0.52 

7deg200vel 0.88 ± 0.04 0.81 ±0.08 0.88 ± 0.04 0.91 ±0.03 5.81 ± 1.11 6.98 ± 1.07 5.90 ± 1.14 7.90 ± 0.38 

8deg200vel 0.91 ± 0.03 0.85 ± 0.06 0.90 ± 0.03 0.93 ± 0.03 6.61 ± 0.76 7.25 ± 0.78 6.59 ±0.81 8.34 ± 0.25 

5deg220vel 0.89 ± 0.04 0.81 ±0.11 0.89 ± 0.04 0.93 ± 0.03 6.02 ± 0.92 6.86 ± 0.88 6.08 ± 0.94 7.24 ± 0.93 

6deg220vel 0.91 ± 0.03 0.81 ±0.08 0.91 ± 0.03 0.93 ± 0.02 5.89 ± 1.02 6.76 ± 1.00 5.92 ± 1.02 7.17 ±1.06 

7deg220vel 0.93 ± 0.02 0.85 ± 0.06 0.93 ± 0.02 0.92 ± 0.04 6.56 ± 0.77 7.20 ± 0.86 6.58 ± 0.80 7.61 ± 0.87 

8deg220vel 0.94 ± 0.02 0.88 ± 0.05 0.94 ± 0.01 0.93 ± 0.02 7.08 ± 0.64 7.95 ± 0.64 7.09 ± 0.66 8.01 ± 0.80 



5deg240vel 


0.96 


± 0.02 


0.88 ± 


0.09 


0.96 


± 0.02 


0.92 


± 0.04 


6.87 


± 0.51 


7.48 ± 0.42 


6.87 


±0.51 


7.46 


±0.43 


6deg240vel 


0.97 


±0.02 


0.84 ± 


0.07 


0.97 


±0.02 


0.89 


±0.06 


6.86 


±0.50 


7.56 ± 0.46 


6.86 


±0.50 


7.54 


± 0.51 


7deg240vel 


0.98 


±0.01 


0.85 ± 


0.06 


0.98 


±0.01 


0.91 


±0.04 


7.46 


±0.25 


7.75 ± 0.39 


7.46 


±0.25 


7.79 


±0.39 


8deg240vel 


0.97 


±0.02 


0.90 ± 


0.04 


0.98 


±0.01 


0.92 


±0.03 


7.58 


±0.21 


7.94 ± 0.39 


7.58 


±0.24 


7.85 


±0.40 



e: mean orbit eccentricity, uncertainties express the standard deviations of the distributions. 

D: mean initial distances in kpc from the infalling galaxy's centre of density, uncertainties as before. 
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Table B.3. Equal-mass merger, infalling galaxy retrograde: particle numbers. Parameters defined as with Table IB. fl 



Model 


^'ini [VparabJ 


rmin [kpc] 


Steps analysed 




N* m 


pall 
retro 


pfar 
retro 


^y60-250kpc 


fin — is i ■„.. — 

pOU— ZjUKpc 

retro 


0deg050vel 


0.5 


0.08 


15 


10355 


47 


0.99 


0.96 


61 


0.48 


2.5deg050vel 


0.5 


1.00 


15 


16314 


7 


0.97 


0.57 


99 


0.68 


5deg050vel 


0.5 


1.93 


14 


19276 


10 


0.93 


0.40 


178 


0.60 


7.5deg050vel 


0.5 


2.80 


15 


26740 


41 


0.93 


0.51 


123 


0.56 


10deg050vel 


0.5 


3.67 


15 


25168 


326 


0.92 


0.79 


195 


0.49 


15deg050vel 


0.5 


5.70 


15 


17146 


1094 


0.96 


0.67 


545 


0.61 


20deg050vel 


0.5 


8.09 


14 


15572 


3842 


0.95 


0.90 


2959 


0.89 


OdeglOOvel 


1.0 


0.16 


15 


10705 


770 


0.88 


0.82 


1973 


0.57 


2.5degl00vel 


1.0 


1.84 


15 


13672 


2657 


0.88 


0.73 


3079 


0.81 


5degl00vel 


1.0 


3.80 


15 


12980 


4313 


0.78 


0.55 


4271 


0.54 


7.5degl00vel 


1.0 


5.97 


15 


15574 


6050 


0.94 


0.91 


4322 


0.76 


1 OdeglOOvel 


1.0 


8.64 


15 


12442 


3780 


0.99 


1.00 


1704 


0.92 


15degl00vel 


1.0 


15.60 


13 


5848 


188 


0.63 


0.40 


516 


0.23 


20degl00vel 


1.0 


25.06 


8 out of 8 


13687 


579 


0.83 


0.63 


267 


0.43 


0degl50vel 


1.5 


0.12 


14 


11907 


3244 


0.88 


0.95 


2668 


0.79 


2.5degl50vel 


1.5 


2.78 


14 


17636 


4061 


0.95 


0.96 


2991 


0.86 


5degl50vel 


1.5 


6.01 


13 


20405 


2932 


0.98 


0.99 


2214 


0.88 


7.5degl50vel 


1.5 


10.15 


9 out of 9 


9008 


1346 


0.96 


1.00 


2299 


0.72 



Table B.4. Equal-mass merger, infalling galaxy retrograde: properties of the two populations. Parameters defined as with Table 



Model gr gr e ^ ^ gar [kpc] gar [kpc] gg [kpc] ^ [kpc] 

0deg050vel 0.78 ± 0.08 0.46 ± 0.09 0.81 ± 0.04 0.84 ± 0.05 7.36 ± 0.75 5.87 ± 0.97 8.67 ± 0.01 6.65 ± 0.56 

2.5deg050vel 0.76 ± 0.08 0.54 ± 0.14 0.86 ± 0.04 0.79 ± 0.08 7.07 ± 0.89 5.78 ± 1.03 7.17 ± 2.12 6.85 ± 1.75 

5deg050vel 0.77 ± 0.08 0.56 ±0.13 0.88 ± 0.04 0.80 ±0.12 7.18 ± 0.72 5.82 ± 1.08 7.55 ± 1.94 7.09 ± 2.00 

7.5deg050vel 0.79 ± 0.09 0.64 ±0.14 0.87 ± 0.04 0.83 ± 0.04 7.25 ± 0.86 6.05 ± 1.18 7.29 ± 1.10 7.24 ±0.91 

10deg050vel 0.78 ± 0.09 0.67 ±0.13 0.85 ± 0.04 0.85 ± 0.05 7.23 ± 0.74 6.46 ±1.17 7.13 ±0.81 6.99 ± 0.86 

15deg050vel 0.83 ±0.13 0.70 ±0.10 0.92 ± 0.05 0.83 ± 0.04 6.78 ± 0.93 6.83 ± 1.08 6.60 ± 0.77 7.18 ±0.57 

20deg050vel 0.88 ±0.14 0.86 ± 0.07 0.95 ± 0.06 0.92 ± 0.04 6.81 ± 1.06 6.33 ± 1.07 6.78 ± 0.98 6.79 ± 0.73 



OdeglOOvel 


0.74 = 


t 0.11 


0.66 = 


t 0.15 


0.88 = 


tO.04 


0.79 = 


= 0.07 


6.80 = 


= 0.86 


6.74 =1 


= 1.20 


7.17 = 


t 0.97 


6.56 ± 


1.19 


2.5degl00vel 


0.78 = 


t0.08 


0.65 = 


t 0.13 


0.83 = 


t 0.05 


0.82 = 


tO.04 


6.97 d 


= 0.73 


6.59=1 


= 1.14 


7.05 = 


t 0.71 


5.98 ± 


1.00 


5degl00vel 


0.81 = 


t0.12 


0.71 = 


t0.09 


0.86 = 


t 0.06 


0.82 = 


tO.04 


6.94=1 


= 0.89 


6.43 d 


= 1.04 


6.97 = 


t0.87 


6.42 ± 


0.87 


7.5degl00vel 


0.80 = 


h0.16 


0.78 = 


t0.09 


0.90 = 


hO.08 


0.84 = 


t0.05 


7.17 d 


= 0.65 


5.97 d 


= 1.09 


7.07 = 


t 0.67 


6.33 ± 


1.05 


1 OdeglOOvel 


0.86 = 


t0.06 


0.80 = 


t0.09 


0.89 = 


t0.02 


0.87 = 


t0.04 


7.28 d 


= 1.12 


6.42 d 


= 1.03 


5.43 = 


h0.23 


6.49 ± 


1.05 


15degl00vel 


0.82 = 


t 0.08 


0.84 = 


t0.10 


0.95 = 


t0.04 


0.96 = 


t0.02 


5.03 d 


= 0.97 


6.35 d 


= 1.26 


6.55 = 


t 1.18 


5.96 ± 


1.14 


20degl00vel 


0.92 = 


t0.04 


0.82 = 


t0.09 


0.93 = 


hO.04 


0.96 = 


t 0.03 


5.66 = 


= 1.14 


6.84 d 


= 0.98 


6.25 = 


t 1.17 


5.98 ± 


1.06 


0degl50vel 


0.62 = 


t 0.22 


0.73 = 


t 0.12 


0.84 = 


tO.06 


0.76 = 


t 0.09 


7.59 d 


= 0.80 


6.30 d 


= 1.16 


7.69 = 


t 1.11 


6.63 ± 


1.18 


2.5degl50vel 


0.67 = 


t 0.20 


0.73 = 


t 0.11 


0.83 = 


tO.07 


0.78 = 


t 0.06 


7.55 =i 


= 0.88 


6.41 d 


= 1.10 


7.73 = 


t 1.10 


6.64 ± 


0.93 


5degl50vel 


0.76 = 


b0.16 


0.70 = 


t0.14 


0.88 = 


t 0.06 


0.75 = 


t0.10 


7.49 d 


= 0.81 


6.04 d 


= 1.20 


8.16 = 


= 0.46 


6.51 ± 


1.25 


7.5degl50vel 


0.84 = 


t0.09 


0.77 = 


tO.10 


0.93 = 


hO.OO 


0.88 = 


t0.07 


5.00 =: 


= 1.12 


6.26 d 


= 1.14 


6.46 = 


t 1.19 


6.23 ± 


1.23 
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Table B.5. Equal-mass merger, infalling galaxy prograde: particle numbers. Parameters defined as with Table IB. fl 



Model 


v ini [ v parab] 


r mm [kpc] 


Steps analysed 




jyfar 


pall 
retro 


pfar 
retro 


^y60-250kpc 


Ail 1? L, v . 

retro 


0deg050vel 


0.5 


0.06 


15 


12092 


47 


0.02 


0.04 


67 


0.33 


2.5deg050vel 


0.5 


1.01 


15 


10660 


689 


0.05 


0.00 


222 


0.04 


5deg050vel 


0.5 


1.95 


14 


14060 


2144 


0.07 


0.00 


1030 


0.01 


7.5deg050vel 


0.5 


2.90 


15 


18794 


3854 


0.06 


0.03 


2063 


0.01 


10deg050vel 


0.5 


3.85 


15 


22421 


4936 


0.02 


0.02 


3329 


0.01 


15deg050vel 


0.5 


5.91 


15 


37049 


6678 


0.04 


0.02 


4832 


0.05 


20deg050vel 


0.5 


8.12 


14 


38768 


10862 


0.03 


0.00 


5499 


0.01 


OdeglOOvel 


1.0 


0.13 


15 


10642 


690 


0.12 


0.20 


1906 


0.44 


2.5degl00vel 


1.0 


2.00 


15 


4644 


925 


0.21 


0.62 


1231 


0.57 


5degl00vel 


1.0 


3.94 


15 


26051 


9347 


0.09 


0.19 


4440 


0.28 


7.5degl00vel 


1.0 


6.09 


15 


42760 


17322 


0.14 


0.13 


10387 


0.20 


1 OdeglOOvel 


1.0 


8.72 


15 


33369 


18475 


0.04 


0.01 


12315 


0.06 


15degl00vel 


1.0 


15.68 


13 


36641 


6143 


0.02 


0.01 


4415 


0.00 


20degl00vel 


1.0 


25.09 


8 out of 8 


35454 


4620 


0.00 


0.00 


9607 


0.00 


0degl50vel 


1.5 


0.07 


14 


11498 


3363 


0.12 


0.05 


2819 


0.21 


2.5degl50vel 


1.5 


2.91 


14 


12372 


2246 


0.27 


0.25 


3466 


0.42 


5degl50vel 


1.5 


6.14 


13 


26998 


7902 


0.20 


0.23 


8727 


0.39 


7.5degl50vel 


1.5 


10.24 


9 out of 9 


25122 


6423 


0.16 


0.11 


11465 


0.15 



Table B.6. Equal-mass merger, infalling galaxy prograde: properties of the two populations. Parameters defined as with Table 
IB~2l 



Model 


al 
e pr 


[ 


all 
c ret 


ro 


ta 
e pr 


r 


.tar 
retro 


pro 


[kpc] 


-^retri 


, [kpc] 


JJtU 
pro 


[kpc] 


^ retro 


[kpc] 


0deg050vel 


0.48 ± 


0.09 


0.79 ± 


0.08 


0.83 ± 


0.04 


0.82 = 


= 0.04 


5.77 =: 


= 0.98 


7.26 


± 0.71 


6.54 = 


t 0.52 


8.67 = 


tO.01 


2.5deg050vel 


0.52 ± 


0.14 


0.80 ± 


0.07 


0.82 ± 


0.03 






5.79 = 


= 1.01 


6.97 


±0.94 


6.76 = 


t 0.68 






5deg050vel 


0.66 ± 


0.17 


0.84 ± 


0.07 


0.80 ± 


0.04 


0.85 = 


t0.04 


6.34 d 


= 1.13 


5.91 


± 1.12 


6.96 = 


t 0.78 


6.16 = 


t 0.44 


7.5deg050vel 


0.73 ± 


0.13 


0.85 ± 


0.06 


0.76 ± 


0.05 


0.88 = 


= 0.05 


6.71 d 


= 1.00 


5.59 


± 1.04 


7.07 = 


t 0.81 


5.66 = 


t 1.10 


10deg050vel 


0.74 ± 


0.10 


0.87 ± 


0.07 


0.71 ± 


0.05 


0.91 = 


t0.05 


6.66 4 


= 1.02 


5.65 


± 1.07 


7.08 = 


= 0.81 


5.76 = 


t 1.08 


15deg050vel 


0.73 ± 


0.11 


0.79 ± 


0.10 


0.65 ± 


0.07 


0.95 = 


t0.02 


5.89 = 


= 1.32 


6.59 


± 1.11 


7.09 = 


t0.75 


5.19 = 


= 0.78 


20deg050vel 


0.64 ± 


0.16 


0.79 ± 


0.11 


0.67 ± 


0.16 


0.96 = 


= 0.00 


5.86 = 


= 1.21 


6.96 


±0.65 


6.18 = 


t 1.08 


7.81 = 


tO.OO 


OdeglOOvel 


0.66 ± 


0.15 


0.74 ± 


0.11 


0.78 ± 


0.06 


0.88 = 


= 0.04 


6.74 ± 1.19 


6.79 


±0.86 


6.50 = 


t 1.22 


7.21 = 


t 1.00 


2.5degl00vel 


0.66 ± 


0.16 


0.82 ± 


0.10 


0.86 ± 


0.05 


0.88 = 


= 0.04 


6.92 ± 1.28 


6.57 


±0.79 


7.36 = 


= 0.93 


6.58 ± 0.74 


5degl00vel 


0.79 ± 


0.11 


0.89 ± 


0.06 


0.83 ± 


0.05 


0.91 = 


= 0.04 


5.98 = 


= 1.25 


5.99 


± 1.09 


6.88 = 


t 0.71 


6.00 ± 1.10 


7.5degl00vel 


0.76 ± 


0.09 


0.76 ± 


0.18 


0.77 ± 


0.08 


0.94 = 


t0.02 


5.89 d 


= 1.33 


6.02 


± 1.27 


6.75 = 


= 0.95 


5.38 ± 0.86 


1 OdeglOOvel 


0.74 + 


0.10 


0.85 ± 


0.10 


0.74 ± 


0.11 


0.95 = 


= 0.05 


5.98 = 


= 1.34 


5.87 


± 1.36 


6.52 = 


t 1.23 


5.48 = 


t 1.45 


15degl00vel 


0.54 + 


0.14 


0.84 ± 


0.08 


0.67 ± 


0.11 


0.91 = 


= 0.03 


5.76 = 


= 1.22 


7.51 


±0.60 


6.44 = 


t 1.16 


7.93 = 


= 0.31 


20degl00vel 


0.51 ± 


0.17 


0.80 ± 


0.19 


0.56 ± 


0.13 






5.62 = 


= 1.32 


7.01 


± 0.70 


6.21 = 


t 1.22 






0degl50vel 


0.73 ± 


0.12 


0.63 ± 


0.22 


0.75 ± 


0.09 


0.84 = 


t 0.06 


6.31 ± 1.15 


7.60 


± 0.77 


6.61 = 


t 1.18 


7.88 = 


= 0.88 


2.5degl50vel 


0.77 ± 


0.12 


0.67 ± 


0.18 


0.85 ± 


0.06 


0.85 = 


= 0.07 


6.05 ± 1.26 


7.20 


±0.85 


6.44 = 


t 1.20 


7.38 = 


= 0.64 


5degl50vel 


0.74 ± 


0.15 


0.68 ± 


0.19 


0.86 ± 


0.09 


0.79 = 


t 0.13 


5.61 d 


= 1.34 


6.66 


± 1.22 


6.62 = 


t 1.29 


6.90 ± 1.08 


7.5degl50vel 


0.67 ± 


0.17 


0.68 ± 


0.17 


0.75 ± 


0.12 


0.80 = 


t0.ll 


5.68 = 


= 1.22 


6.75 


± 1.35 


5.82 = 


t 1.11 


7.11 = 


t 1.19 



Table B.7. 4-to-l-mass fly-by: particle numbers. Parameters defined as with Table IB. ll 



Model 


^ini [^'parab] 


r mm [kpc] 


Steps analysed 




yyfar 


pall 

retro 


pfar 
retro 


^y60-250kpc 


p60-25()kpc 
retro 


4degl75vel 


1.75 


6.89 


16 


11767 


3813 


0.33 


0.00 


2216 


0.02 


6degl75vel 


1.75 


11.35 


16 


4116 


2026 


0.02 


0.00 


1581 


0.04 


8degl75vel 


1.75 


16.86 


16 


3951 


1746 


0.00 


0.00 


1179 


0.00 


2deg200vel 


2.00 


3.67 


16 


1263 


628 


0.03 


0.05 


458 


0.09 


4deg200vel 


2.00 


7.84 


15 


12689 


4745 


0.50 


0.01 


3203 


0.01 


6deg200vel 


2.00 


13.09 


16 


3295 


2279 


0.08 


0.00 


1660 


0.01 


4deg225vel 


2.25 


8.84 


16 


13942 


6095 


0.64 


0.21 


4677 


0.08 


6deg225vel 


2.25 


14.76 


16 


3684 


2919 


0.04 


0.01 


2303 


0.01 


4deg250vel 


2.50 


9.73 


16 


9266 


5837 


0.94 


0.90 


5049 


0.56 


6deg250vel 


2.50 


16.21 


15 (steps 4k to 18k) 


1265 


1138 


0.17 


0.19 


1674 


0.10 



Pawlowski, Kroupa & de Boer: Making Counter-Orbiting Tidal Debris 
Table B.8. 4-to-l-mass fly-by: properties of the two populations. Parameters defined as with Table |B.2l 



Model 


all 

''pro 


all 

retro 


tar 
pro 


e 


Far 
retro 




[kpc] 


B%L [kpc] 




[kpc] 


B&o [ k P c ] 


4degl75vel 


0.83 


± 0.06 


0.81 ± 


0.06 


0.83 ± 


0.06 


0.88 


± 0.04 


4.06 


± 0.85 


4.68 


± 0.82 


3.98 


± 0.86 


5.90 


±0.48 


6degl75vel 


0.88 


±0.05 


0.84 ± 


0.06 


0.88 ± 


0.04 


0.98 


± 0.00 


4.07 


± 0.81 


5.43 


±0.95 


4.04 


± 0.82 


6.17 


±0.00 


8degl75vel 


0.88 


± 0.05 


0.86 ± 


0.05 


0.91 ± 


0.03 






4.59 


± 0.77 


5.11 


±0.97 


4.41 


± 0.74 






2deg200vel 


0.83 


±0.07 


0.95 ± 


0.02 


0.87 ± 


0.05 


0.95 


± 0.02 


5.68 


± 0.47 


5.69 


± 0.57 


5.69 


± 0.57 


5.67 


±0.57 


4deg200vel 


0.88 


±0.04 


0.82 ± 


0.06 


0.88 ± 


0.04 


0.90 


±0.04 


4.09 


± 0.77 


4.88 


±0.89 


4.07 


±0.78 


5.69 


±0.90 


6deg200vel 


0.93 


±0.04 


0.86 ± 


0.06 


0.94 ± 


0.02 


0.94 


±0.02 


4.43 


± 0.70 


5.32 


±0.68 


4.40 


±0.69 


5.71 


±0.84 


4deg225vel 


0.94 


±0.02 


0.85 ± 


0.06 


0.94 ± 


0.02 


0.91 


± 0.04 


4.26 


±0.68 


4.64 


± 0.74 


4.26 


± 0.67 


4.63 


±0.81 


6deg225vel 


0.94 


±0.03 


0.89 ± 


0.04 


0.95 ± 


0.02 


0.92 


±0.03 


4.48 


±0.74 


6.26 


±0.40 


4.39 


± 0.75 


5.96 


±0.58 


4deg250vel 


0.97 


± 0.01 


0.87 ± 


0.06 


0.97 ± 


0.01 


0.91 


± 0.04 


4.40 


±0.67 


4.98 


±0.67 


4.40 


± 0.67 


4.84 


± 0.74 


6deg250vel 


0.95 


±0.02 


0.94 ± 


0.03 


0.96 ± 


0.01 


0.95 


±0.03 


5.44 


±0.30 


5.87 


±0.38 


5.42 


± 0.31 


5.86 


±0.38 



Table B.9. 4-to-l-mass merger, retrograde infalling galaxy: particle numbers. Parameters defined as with Table IB. ll 



Model 


v ini [ v parab] 


r mm [kpc] 


Steps analysed 


N M 


2yfar 


pall 
retro 


pfar 
retro 


jy60-250kpc 


„60-2J0kpc 
retro 


2.5deg050vel 


0.5 


1.16 


13 


62165 


19714 


0.63 


0.51 


17625 


0.37 


5deg050vel 


0.5 


2.33 


13 


49654 


23292 


0.62 


0.62 


19567 


0.45 


2.5degl00vel 


1.0 


2.28 


13 


62118 


27878 


0.45 


0.58 


28341 


0.53 


5degl00vel 


1.0 


4.76 


13 


28922 


18078 


0.44 


0.48 


20009 


0.38 


7.5degl00vel 


1.0 


7.50 


13 


13263 


8192 


0.23 


0.03 


13521 


0.14 


lOdeglOOvel 


1.0 


11.02 


11 out of 11 


15762 


2067 


0.30 


0.06 


5071 


0.48 



Table B.10. 4-to-l-mass merger, retrograde infalling galaxy: properties of the two populations. Parameters defined as with Table 
lB~2l 



Model 


all 




all 
retro 


e 


tar 
pro 


retro 


D aU 

pro 


[kpc] 


D m 

-^retro 


[kpc] 


jytu 

pro 


[kpc] 


^ retro 


[kpc] 


2.5deg050vel 


0.86 ± 


0.07 


0.86 = 


t 0.06 


0.89 


±0.06 


0.85 = 


t 0.05 


4.08 : 


± 1.03 


3.91 = 


= 0.93 


4.53 : 


±0.89 


4.34 = 


= 0.87 


5deg050vel 


0.83 ± 


0.10 


0.87 = 


t 0.06 


0.88 


± 0.07 


0.87 = 


t 0.06 


4.35 : 


± 1.02 


4.16 = 


= 0.93 


4.63 : 


±0.92 


4.36 = 


= 0.83 


2.5degl00vel 


0.85 ± 


0.08 


0.87 = 


t 0.05 


0.83 


± 0.10 


0.88 = 


t 0.05 


3.46: 


± 1.15 


4.20 = 


= 1.00 


4.40: 


±0.96 


4.55 = 


t 0.91 


5degl00vel 


0.86 ± 


0.09 


0.88 = 


t 0.05 


0.89 


± 0.08 


0.90 = 


t0.04 


4.20: 


± 1.01 


4.62 =i 


= 0.94 


4.53 : 


±0.93 


4.96 = 


t0.73 


7.5degl00vel 


0.87 ± 


0.12 


0.86 = 


t0.05 


0.92 


±0.05 


0.92 = 


t0.03 


4.35 : 


± 1.05 


4.26=i 


= 0.86 


4.42 : 


± 1.01 


5.55 = 


t0.50 


lOdeglOOvel 


0.76 ± 


0.11 


0.84 = 


t0.04 


0.90 


±0.04 


0.90 = 


t0.03 


4.25 : 


± 1.04 


4.26 = 


= 0.86 


4.53 : 


±0.88 


5.10 = 


t0.99 



Table B.ll. 4-to-l-mass merger, prograde infalling galaxy: particle numbers. Parameters defined as with Table IB. ll 



Model 


Vini [Vparab] 


''mill [kpc] 


Steps analysed 


N M 


rV far 


pall 
retro 


nfar 
retro 


jy60-250kpc 


„6()-2J0kpc 
'retro 


2.5deg050vel 


0.5 


1.29 


13 


42268 


10307 


0.44 


0.66 


14741 


0.66 


5deg050vel 


0.5 


2.44 


13 


47474 


12688 


0.26 


0.39 


16696 


0.50 


2.5degl00vel 


1.0 


2.38 


13 


70820 


26464 


0.28 


0.53 


27727 


0.61 


5degl00vel 


1.0 


4.79 


13 


62362 


31521 


0.06 


0.09 


29706 


0.10 


7.5degl00vel 


1.0 


7.49 


12 


46993 


29979 


0.01 


0.00 


28675 


0.04 


lOdeglOOvel 


1.0 


11.09 


11 out of 11 


33946 


19700 


0.05 


0.01 


23483 


0.05 



Table B.12. 4-to-l-mass merger, prograde infalling galaxy: properties of the two populations. Parameters defined as with Table 
IB~2l 



Model 


all 
°pro 


all 
c retro 


tar 
^pro 


tar 
retro 


D aU 

^pro 


[kpc] 


D all 

-^retro 


[kpc] 


D m 

pro 


[kpc] 


^ retro 


[kpc] 


2.5deg050vel 


0.82 =i 


= 0.11 


0.89 = 


t 0.07 


0.86 = 


= 0.09 


0.94 = 


= 0.04 


3.61 


± 1.05 


3.92 = 


= 0.83 


4.52 : 


± 1.05 


4.18 = 


t 0.77 


5deg050vel 


0.86 = 


= 0.08 


0.90 = 


t 0.06 


0.90 = 


= 0.06 


0.93 = 


= 0.04 


3.46 


±0.99 


3.96 = 


t 0.74 


3.92: 


± 1.16 


4.12 = 


= 0.75 


2.5degl00vel 


0.84 d 


= 0.07 


0.90 = 


t 0.07 


0.90 = 


= 0.06 


0.93 = 


t 0.06 


3.36 


± 1.04 


4.05 = 


t0.85 


3.86: 


± 1.08 


4.04 = 


t 0.81 


5degl00vel 


0.85 d 


= 0.09 


0.89 = 


t 0.09 


0.91 = 


= 0.05 


0.92 = 


t 0.04 


3.35 


±0.92 


4.10 = 


= 0.98 


3.64 


±0.89 


3.98 = 


= 0.89 


7.5degl00vel 


0.84 = 


= 0.07 


0.88 = 


= 0.09 


0.87 d 


= 0.06 


0.95 = 


= 0.03 


3.22 


±0.89 


3.85 = 


t 1.03 


3.37 : 


±0.89 


4.01 = 


t 1.04 


lOdeglOOvel 


0.82 = 


= 0.08 


0.84 = 


tO.09 


0.84 = 


= 0.06 


0.89 = 


= 0.08 


3.41 


±0.93 


3.19 = 


t 1.05 


3.43 : 


±0.86 


3.68 = 


= 0.80 



26 



